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Abstract. We consider some interacting particle processes with long-range 
dynamics: the zero-range and exclusion processes with long jumps. We prove 
that the hydrodynamic limit of these processes corresponds to a (non-linear 
in general) fractional heat equation. The scaling in this case is superdiffusive. 
In addition, we discuss a central limit theorem for a tagged particle on the 
zero-range process and existence and uniqueness of solutions of the Cauchy 
problem for the fractional heat equation. 



1. Introduction 

Interacting particle systems have been the subject of intense study for the last 40 
years. This is due to the fact that, in one hand, they present many of the collective 
features found in real physical systems, and in the other hand they are, up to 
some extent, mathematically tractable. The rigorous study of interacting particle 
systems has lead in many cases to a more detailed understanding of the microscopic 
mechanism behind those collective phenomena. We refer to |^ for further references 
and to [3] for a recent result which we think is a good example of the success this 
plan. 

Since their introduction by Spitzer, the zero-range process and the exclusion 
process have been among the most studied interacting particle systems, and they 
have served as a test field for new mathematical and physical ideas. 

During the last years, and specially due to applications in finance, fractional 
Laplacians L — —(—A)*, s € (0, 1) and their probabilistic counterparts, the Levy 
processes, have received grown attention. Some of the key properties of Levy pro- 
cesses are the presence of jumps, the lack (at least in the cases associated to L) of 
bounded variance and a super-diffusive behavior. 

From the point of view of statistical mechanics, it is desirable to have a deriva- 
tion of the partial differential equations ruling the evolution of these super-diffusive 
systems, which usually involve the fractional Laplacian, from microscopic models. 
We are not aware of any work on this direction, so we have decided to provide 
in this article such a derivation. Our aim is to obtain the hydrodynamic limit of 
interacting particle systems on which particle may perform long jumps, in the con- 
text of the exclusion and zero-range processes. For these models, the corresponding 
hydrodynamic equation is given by a fractional (non linear in general) heat equa- 
tion of the form dtu = Lu, where L is the generator of a symmetric, a-stable Levy 
process and in particular it includes the fractional Laplacians — (—A)'' as particular 
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cases. We believe that this equation wiU emerge as weh as the hydrodynamic hmit 
of other particle systems which are superdifFusive in nature, like heat conduction 
models with conservative noise in dimension d = I, Ornstein-Uhlenbeck particles 
in laminar flow, etc. 

This article is organized as follows. In Section [2] we give detailed definitions 
of the models considered here. In Section [3] we prove our main result for the 
exclusion process with long jumps. At the current state of art, this result is more 
or less elementary. We included a complete proof for two reasons. In general 
terms, the theory of hydrodynamic limits of particle systems is a hard subject for 
the non-expert reader, since it involves a mixture of purely probabilistic ideas and 
ideas coming directly from analysis and PDE theory. Therefore, we believe that 
it is a good way to introduce the subject in a more elementary way, simplifying 
later the exposition for the zero-range process. And secondly, the results of this 
section were used without proof in 6 to obtain an invariance principle for a tagged 
particle in the exclusion process with long jumps. In Section [4] we prove the main 
theorem for the zero-range process, leaving various technical results for the following 
sections. In Section [5] we prove the key technical input, known in the literature as 
the replacement lemma. In Section [6] we prove the so-called moving particle lemma, 
which we believe is the main feature that differentiates super-diffusive systems from 
the most studied diffusive systems, and the new proof of this lemma is the main 
technical novelty of this work. Although a portion of the exposition is by now 
classical, the tools needed have been gathered from many different places, and we 
have decided to include detailed proofs of most of the propositions taken from 
elsewhere in order to keep the exposition as clear as possible. We invite the most 
specialized reader to skip the more standard parts of the exposition. 

In Section [7] we prove the energy estimate, which is crucial in order to obtain a 
uniqueness criterion for the hydrodynamic equation. In this section we introduce a 
variational formula for the Fisher information for the fractional heat equation which 
seems to be new in the literature. This formula involves a natural generalization of 
the space of test functions to antisymmetric functions of two space coordinates. In 
the derivation of hydrodynamic limits of particle systems, a key analytical input is 
an existence result for the corresponding hydrodynamic equation. After consulting 
some experts in the field, it seems that the required uniqueness results are not avail- 
able in the literature, so in Section [8] we obtain some uniqueness results tailored to 
our needs. Those results may be of independent interest, and they are independent 
of the rest of paper. The proof in the linear case is due to Luis Silvestre. 

In Section [9] we prove an invariance principle for a tagged particle in the one- 
dimensional zero-range process, following recent results in [6], [7]. The limiting 
process is a time-inhomogeneous process of independent increments, related to the 
solution of the hydrodynamic equation. In Appendix [A] we explain how to deal 
with the borderline case a — 2, which leads to the usual heat equation, but with 
a superdiffusive time scaling. The hydrodynamic limit of the zero-range process in 
its full generality is proved only under a restrictive attractiveness condition for the 
system. However, we point out that it is only at the level of the uniqueness criterion 
for the hydrodynamic equation that this result is needed. Certainly most of the 
proofs can be obtained in a simpler way using coupling arguments, available only for 
attractive systems. In this article, attractiveness is only used to obtain a bound for 
the energy of the solutions in terms of their Fisher information. We conjecture that. 
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like in the diffusive case, uniqueness under the weaker bound on Fisher information 
holds. In Appendix iBl we explain how to consider unbounded initial profiles in the 
formulation of the hydrodynamic limit for attractive systems. In Appendix [C] we 
explain how to handle general bounded initial profiles. Both results are based on 
coupling tecniques which are more or less standard in the literature. A combination 
of both appendices can be used to handle initial conditions of the form uo{x)+vq{x), 
where uq is bounded and vq is in £^(M''). 



2. The models 

Consider the integer lattice 'Z'^, d > 1 and let p : Z"* [0, oo) be such that 
Szgz<iP(-^) = P* < +00. We call p(-) the transition rate. The exclusion process 
associated to p{-) is the Markov process rjt defined in Q^x — {0, 1}^'' and generated 
by the operator 

x,y£Z'' 

Here 77 denotes a generic element of flex, f ■ ^ex ^ M is a local function, that 
is, it depends on the value of 77(2) for only a finite number of points z G Z'', and 

ri{z), z^x,y. 

For details about the construction and properties of the process rjt, we refer 
to [12]. The dynamics of this process is easy to describe. Initially, particles are 
distributed in Z*^ in such a way that there is at most one particle per site. Each 
particle, independent of the other particles, waits an exponential time of rate p* , at 
the end of which it picks a site y S Z'^ with probability p{y — x)/p* , where x is the 
current position of the particle. If the site y is empty, the particle jumps from x 
to y. Otherwise the particle stays at x. In any case, a new exponential time starts 
afresh. Notice that for initial configurations with a finite number of particles, the 
process rjt is just a system of independent random walks with transition rate p{-), 
conditioned to have at most one particle per site. When the number of particles is 
infinite, the construction of the process has been carried out by Liggett jl 1] . 

Denote by Nq the set of non-negative integers and let g : Nq ^ [0,oo) be a 
function such that g{0) — 0. We assume that g{n) > for n > 0. The zero-range 
process with interaction rate g{-) and transition rate p{-) is the Markov process 
defined in fl^^, = (Nq)^'' and genetared by the operator 

i^v(o= E p(2^-^)5(cw)[/r'^)-/(o]- 

x.yel.'^ 

Here ^ is a generic element of 1^°^, / : il^r ^ K is a "suitable" function and 

{£,{x) — 1, z — X 
e(y) + i, z^y 
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The exact meaning of "suitable" depends on the choice of p(-), but for local 
functions /, a Lipschitz condition of the form 



for any ^, x, y will be sufficient. 

The dynamics of this process is the following: at a site x € Z'' and independently 
for each site, the particles wait an exponential time of rate p*g{£,{x)), where ^(x) 
denotes the number of particles at site x. At the end of this exponential time, one of 
the particles at x jumps to a site y, randomly chosen with probability p{y — x) /p* . 
Then, a new exponential time starts afresh. Notice that the particles interact 
between them only when they are at the same site, explaining the denomination 
"zero-range process". The case g{n) = n corresponds to a system of independent 
random walks in Z"* with transition rate p{-). 

Differently from the exclusion process, the zero-range process can present explo- 
sions if the initial configuration of particles has too many particles at infinity. In 
this article, we assume that the interaction rate satisfies the Lipschitz condition 
stated above: there is a finite constant k > such that \g{n + 1) — g{n)\ < k for 
any n > 0. We refer to this condition as (LG). Under (LG), the process is well- 
defined for any bounded initial condition and also a.s. with respect to any initial 
measure /i in 17°^ such that sup^. £'p[^(x)] < -t-oo. Fore a more detailed discussion 
about this point, the construction of the process and related topics, we refer to [1]. 



2.1. Invariant measures. Let p be a fixed number in [0,1]. Define fip as the 
product measure in fl^x with one-site marginals satisfying 



It is well known that the measures {/ip; p G [0, 1]} are invariant under the 
evolution of rjt. When p{-) is symmetric, the measures {lJ-p}p are also reversible. 
And if the transition rate p(-) is irreducible, then the measures {A*p}p are also 
ergodic. Notice that J r]{x)dfip = p, that is, the density of particles per site is equal 
to p for the measure jjip. The fact that there exists a family {pp\p of invariant 
measures parametrized by the density o particles reflects the fact that the dynamics 
is conservative: particles are neither created nor destroyed by the dynamics. 

For the zero-range process there also exists a family of invariant measures of 
product form [T]. Take > and define v^j, as the measure in 57°,, of marginals 
given by 



where ^(fc)! — g{l) ■ ■ ■ g{k) for fc > 0, 5(0)! — 1 and Z{(j)) is the normalization 
constant 



Notice that Z((j)) is an increasing function of (j). Therefore, there is a maximal 
value (f>c (perhaps equal to +00) such that Z{(j)) is finite for (j) < (j)c and Z{(j)) is 
infinite for (p > (pc- If there is a positive constant eq such that g{n) > eg for any 
n > 0, 0c > and Z{(j)c) = +00. The family of measures {P^; (j) < 0c} is invariant 
under the dynamics of ^f. When the transition probability p{-) is symmetric, the 



sl.l 



fJ-pivi^) = 1) = 1 - Pp{v{x) = 0) = p. 
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measures are reversible. And when p{-) is irreducible, the measures u^f, are 
ergodic. 

The zero-range process also conserves the number of particles, so it would be 
more natural to parametrize the invariant measures by the density of particles per 
site. Define p{4>) = J £^{x)i'p. It is not difficult to see that p{(f>) is a differentiable, 
strictly increasing function from [0, 0c) to [0, oo). Let us write pc = lim^^;^^ Plv)- 
The inverse function 0(p) is well defined for any p £ [0,pc). We define Vp = v^^py 
Notice that 0(p) = J g{S_{x))vp. 

The occupation variables ^{x) and the interaction rates g{£,{x)), x £ TL'^ have 
exponential moments of sufficiently small order. In fact, for e R, 



and we conclude that Mp{9) is finite for 6 < log((/)c/0). From (LG), we obtain the 
bound g{n) < nn, and therefore / e^3^^^^'>'>dP^ is finite if 6* < log(0c/0). We 
say that the interaction rate g{-) satisfies (FEM) if MpiO) is finite for every 6, p> 
0. A simple application of the ratio test shows that a non-decreasing interaction 
rate g{-) satisfies (FEM) if and only if 4>c — +oo, and equivalently if and only if 
lim„^oo = +00 (this last limit always exists, since g{-) is non-decreasing). We 
say that g{-) satisfies (B) if g{-) is non-decreasing and if it does not satisfy (FEM). 
Therefore, g{-) satisfies (B) if and only if g{-) is non-decreasing and bounded. 

We already mentioned that the process is well defined a.s. with respect to an 
initial measure p such that sup^ i?^[^(a;)] < -l-oo. Here we quote a more precise 
statement, due to Andjel 1 : 

Proposition 2.1. For any transition rate p{-), there exists a function ct : Z'' — > 
[0, oo) such that: 

i) EzCr(z) < -l-oo, 

ii) The zero-range process is a strong Markov process when defined on the 
set 

r!,, = {e e ^ eW'T(z) < +oo}. 

Notice that flzr has full measure for any of the invariant measures Vp, p > 0. 
From now on, we will always define the process in ^zr- 

Remark 2.2. For x e M'^, we denote by \ \x\\ the Euclidean norm {xl + - ■ --l-x^)^/^ 
of X, and by \x\ the supremum norm sup{|a;i|, . . . , Ixdl}- For a given Polish space 
£, we denote by I?([0,oo),£) ('D{[0,T],£) resp.) the space of cddldg functions 
f : [0,oo) — !■ £ (f : [0,T] ^ £ resp.) equipped with the J -Skorohod topology. 

2.2. Homogeneous transition rates. When the transition rate p(-) has mean 
zero and finite range, that is, when zp{z) = and p[z) = for z big enough, it 
is well known that the hydrodynamic limit of the processes r/j, is diffusive and 
given by the heat equation in the case of the exclusion process rjt and by a nonlinear 
heat equation for the process (see [5] and the references therein). In the case 
zp{z) = TO 7^ 0, the hydrodynamic limit appears in the hyperbolic scale and 
corresponds to a conservation law. Therefore, in order to obtain a superdiffusive 
scaling limit, it is necessary for the transition rate p{-) to have arbitrarily long 
jumps. Probably, the most natural choice for such a p{-) should be to take p{x) = 
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where ||a;|| — (xf + • • • + x^^^l"^ is the Eudidean norm of x and a > is 
an arbitrary constant. Notice that p* < +00 due to the condition a > 0. Now we 
describe a broader class of transition rates for which l/||a;|| will be a particular 
case. A function h {0} — > R is said to be homogeneous of degree e M if 

for any x £ M.'^ \ {0} and any A > 0, h{Xx) = X^h{x). We say that a transition 
rate p(-) is homogeneous, regular of degree /? if there is a function h : {0} ^ E 
homogeneous of degree continuous and strictly positive such that p{x) = h{x) 
for any a; G Z'^ \ {0}. 

We assume the following hypothesis on p(-): 

(P): The transition rate p{-) is symmetric and homogeneous of degree — + 
a), with a € (0,2). 

The restriction a > comes from the fact that '^^p{z) < +cx), and in particular 
the function h{-) has to be integrable outside a ball around the origin. The restric- 
tion a < 2 comes from the fact that for a > 2, z'^p{z) < +00 and in that case 
the hydrodynamic limits of rjt and are still diffusive. The boundary case a = 2 
is special, since the hydrodynamic limit is expected to be the usual heat equation, 
but the scaling contains a logarithmic correction (see Appendix |^ . 

Notice that in d = 1, the unique transition rates p(-) satisfying (P) are p{z) — 
c/||z||^+", c > 0. In fact, the class of transition rates p{-) satisfying (P) is home- 
omorphic to the class of measures m in the sphere S''^^ = {x G R'*; ||a;|| = 1} 
with a continuous, strictly positive density with respect to Lebesgue measure and 
satisfying mi A) = m{—A) for any Borel set A G We call the processes rjt 

and associated to a transition rate p(-) satisfying (P) the exclusion process and 
zero-range process with long jumps. 

2.3. The hydrodynamic limit for rjf Let p{-) be a transition rate satisfying 
(P) and let /i(-), —{d + a) be the corresponding homogeneous function and degree. 
These parameters will be fixed throughout the rest of this article. We start defining 
the pseudo-differential operator 

^Fix) = i / h{y){Fix + y) + Fix - y) - 2F{x))dy 

for F E C^(M'*), the set of twice-continuously differentiable functions of bounded 
support. Since the function F is bounded, the integral is absolutely convergent 
outside a ball around the origin (at this point we need a > 0). And using a second- 
order Taylor expansion of F around x we see that the integral is also absolutely 
convergent around (here we need a < 2). In fact, CF : R*^ ^ M'* is a continuous 
function. Notice that in the case h{y) — c/||y||''+", the operator £ corresponds to 
a constant multiple of the fractional Laplacian — (— A)"/^. 

Let Up : K'' ^ R be a measurable function. We say that a sequence of probability 
measures in fJea; is associated to the initial profile uq if for any function F 

in the set Cc(R'') of continuous functions of bounded support and every e > 0, 

lim ^"|?7; n^'' ^ ri{z)F{z/n) — / uo{x)F{x)dx > e| 0. 

We adopt the same definition for a sequence of measures {i^"}n in ^zr^ exchang- 
ing ?7 by ^ in the previous relation. Notice that in order to have a sequence of mea- 
sures {/x"}„ associated to an initial profile uo, it is necessary to have < uo{x) < 1 
for any a; G R'^. 
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Fix r > 0. We say that a measurable function u : [0, T] x M'^ — > M is a weak 
solution of the Cauchy problem 

J dtut = Cut , . 

if llulloo.T = sup{M(i, a;); t G [0, r],a; e W^} is finite and for any smooth function 
G : [0, oo) X K'^ ^ R of compact support we have 

/ u{T,x)G{T,x)dx— / uq{x)G{0, x)dx~ / / 

{t,x){dt+C}G{t,x)dxdt = 0. 

At this point we are ready to define what we mean by the hydrodynamic limit 
of the process 77^. 

I tl I Theorem 2.3. Let uq '■ M'' — > [0, 1] be an initial profile, and let {/i"} be associated 
to uq. Define 77" as the speeded-up process rjtn" starting from the initial measure 
/i" and let be the distribution in flex of r/" . Then, the sequence {A*"(i)}n is 

associated to the function u(t,-), where u(t,x) is the unique weak solution of the 
Cauchy problem 12.1]) . 

Equation (|2.ip is known as the hydrodynamic limit of the process ry". It is 
important to notice the superdiffusive scaling n" in this theorem. For a given 
measurable initial profile uq : M'' ^ [0, 1], it is not difficult to construct a sequence 
of measures {/i"}™ associated to it. In fact, it is enough to consider the product 
measures fi"" defined by 

m"{77(z) = 1} = 1-^"{77(z) = 0}=<(z)=: J uo{x)dx. (2.2) 

Let A^ + (R'^) be the space of non-negative, Radon measures in R''. Let the 
empirical density be the random measure in A^ + (R'') defined by 

7r"(c?x) = n^"^ ^ T]{z)S^/n{dx), 

where Sx is the Dirac-5 distribution at x. It is not difficult to see that the sequence 
is associated to uo if and only if the sequence {7r"(da;)}„ of random measures 
satisfies a weak law of large numbers with respect to the weak topology in A4+(R''), 
with limit measure uo{x)dx. 

1.4 I 

2.4. Atractiveness and stochastic domination. Let ^, ^' be two elements of 
^Izr- We say that ^ ^ ^' if and only if ^{x) < ^'(x) for any x £ 71^ . The relation < 
defines a partial order in ^^.r ■ 

We say the the process r\t is attractive if for any two initial configurations of 
particles ^ in with ^ there exists a process {(,l,£,t) in ^zr x ^zr such 
that the evolution of i = 1,2 corresponds to a zero-range process starting from 

and ^ a.s. for any t>0. 

It is not hard to see that the process is atractive when the interaction rate g{-) 
is non-decreasing. The construction of the process {£,l,£,t) is simple. First we take 
a zero-range process ^} starting from Fix now a realization of the process and 
let 6t be the zero-range process with initial configuration S defined by S{x) = ^^(a;) — 
^^{x) and inhomogeneous interaction rate gx,t{5) = g{6{x)+£,}{x)) — g{£,}{x)). The 
process St is well defined exactly due to the fact that g{-) is non-decreasing. Then 
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we define by taking ^i{x) = £,t{x) + 5t{x). For a justification of the validity of 
tliis construction, see Section 4 of [1]. 

Let v' be two probability measures in ^zr- We say that v is stochastically 
dominated by v' and we write v ^ v' , if there exists a measure {> in x flzr such 
that = j>(-, ri^r), t^'(-) i>(^^^r, ■) and ^ ^ i 

We say that a function F : ^ M is non-decreasing if F(^) < whenever 
C ^ C- Taking the expectation of F{^^) — with respect to i>, it is easy to see 

that if ^ t^', then 

J Fdv < J Fdv' 

for any non-decreasing function F. In fact, the validity of this last relation for 
any non-decreasing, bounded function F is sometimes used as the definition of 
stochastic domination. 

The following property of the process S^t follows easily from its atractiveness. 

Proposition 2.4. Let , he two probability measures defined in fl^r with ^ 
v'^ . Denote by v^{t) the distribution at time t of the zero-range process starting 
from the initial measure v\ i = 1,2. Then ^^{t) is stochastically dominated by 
y\t). 

For any < (? , we have Vp\ -< Vpi (see Section 1 of [S]). Here we give a quick 
probabilistic proof of this fact. It is enough to prove that qp\ -< qpi, where qp is the 
probability measure in No given by the one-site marginal of Vp. Notice that qp is the 
unique invariant measure of the ladder process in Nq defined as follows. A particle 
at site n > goes down with exponential rate g{n) and goes up with exponential 
rate (/'(p). Now take two particles X^, evolving in No as follows. At time t = Q 
they both start at n = 0. When they are not together they evolve independently, 
following a ladder process of rate i — 1,2. When they are together, say at 

site n, they wait an exponential time of rate g{n) at the end of which they go one 
step down together. They also wait an exponential time of rate 4>{p^), at the end of 
which particle X"^ goes one step up, and particle X^ goes up together with particle 
X"^ with probability (j){p^)/(j){p'^) and it stays at n with probability 1 — 
It is clear that the law of X^ corresponds to a ladder process of rate (/'(p*), i — 1,2 
and X^{t) < X'^{t) for any t > 0. The process {X^{t), X'^{t)) constructed in this 
way is recurrent. In fact, X'^{t) is recurrent, so the return time to is a.s. finite 
and has finite mean. But each time X'^{t) returns to 0, X^{t) also returns to 0. 
Therefore, the return time of {X'^{t), X'^{t)) to (0,0) is a.s. finite and has finite 
mean. We conclude that {X^{t),X^{t)) is recurrent, and therefore it has a unique 
invariant measure. For a process with only one invariant measure, the Cesaro means 
of the distributions at time t of the process converge to the invariant measure. In 
our case, the Cesaro means of the distributions of {X^ (t) , X'^ (t)) converge to a 
probability measure q with marginals qpi and satisfying q{x^ < a:^) = 1. 

2.5. Hydrodynamic limit for ^f. Before stating our result about the hydrody- 
namic limit of the zero-range process ^t, we need some definitions. Let /x, u be two 
probability measures defined in some measurable space ft. The relative entropy of 
p with respect to is defined by 

J dv dv 

if p is absolutely continuous with respect to v, and H{p\i') — +oo otherwise. 
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Remember the definition of (t>{p) in Section 12.11 For a function M : Krf ^ K we 
define the energy form £{u, u) as 

^ {u, u) = / / h{y - x) {u{y) ~ u{x)) ^dxdy. 



2 . . 

For any two functions u, v such that £{u,u) < +oo, £{v,v) < +oo, we define 

'^("'«) = ^ J J Hy - x){u{y) ~ u{x)){v{y) - v{x))dxdy. 

Notice that for functions F, G in C^{R'^), £{F,G) = - / F{x)CG{x)dx. Fix a 
reference density p > 0. Define the entropy Ti. : [0,oo) — > [0,oo) as 

H{a) = j\og{^^)dx. (2.3) 

Fix T > and let : R'' R be a bounded function satisfying J \uq{x) — p\dx < 
+00. We say that a measurable function u : R.'^ x [0,T] ^ M is an energy solution 
of the Cauchy problem 

rn*^ = '^'^^"^ (2.4) [7^ 

■u(0, •) = Mo(-) I 1 

if: 

i) ut(-) = u(',i) has finite entropy, that is, 

y J H{u{x,t))dxdt < +00, (2.5) | ec2.4 

ii) ut satisfies the energy estimate 

rT 

£{(t){ut), (j){ut))dt < +00, (2.6) I ec2.3 



iii) for any smooth function G : K.'' x [0, T] ^ M with a compact support 
contained in R'^ x [0,r), 

(t>{u{x,t))CG{x,t) +u{x,t)dtG{x,t)^dxdt + / G{x,0)uo{x)dx ^ 0. (2.7) 



Theorem 2.5. Let Wq : K"* ^ [0,oo) 6e a measurable, bounded initial profile with 
j'{uQ{x) — u)'^dx < +00 for some constant u > and let {j^"}n be a sequence 
of probability measures in fl^r associated to uq. Let ^" be the zero-range process 
S^tn'^ starting from i^" and let v"{t) be the distribution in flzr of . Assume that 
the interaction rate g{-) is non- decreasing and that there exist positive and finite 
constants p, K such that 

(H) For any n>0, H{v'''\vp) < Kn'^ . 

If the interaction rate g(-) satisfies (FEM), also assume that there is a constant 
p' > p such that 

(C) The measures i^" are stochastically dominated by Vpi . 

Then {i^"(t)}„ is associated to the function u{-,t), where u(x,t) is the unique 
energy solution of ^2.4\ )- 



ec2.2 
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We say that equation ()2.4p is the hydrodynamic limit of the zero-range process ■ 
Remember that the function 4'{p) appearing in the hydrodynamic equation is equal 
to the expectation of the interaction rate g(-) with respect to the invariant measure 
Vp. For an integrable initial profile uq, bounded in the complement of some ball 
around the origin, it is not difficult to see that the product measures defined 
in flzr by substituting the Bernoulli marginals in (|2.2p by the measures qu^(z) 
are associated to uq{-). These measures also satisfy hypothesis (H) if / \uo{x) — 
p\dx < +00. When (H) is satisfied, we necessarily have u = p. Hypothesis (C) 
is quite restrictive: it is satisfied by the measures v"' defined above if and only 
if sup^Mo(a;) < p' . In Appendix IB] we explain how to get rid of hypothesis (C). 
Perhaps the most restrictive assumption is being non-decreasing, since this is an 
assumption on the dynamics and not on the initial profile. We will see that the only 
point where we need g{-) to be non-decreasing is to obtain enough conditions on 
the limiting profile u{x,t) in order to guarantee uniqueness of the hydrodynamic 
equation (|2.4p . In Appendix iBl we state a result which does not require g{-) to 
be non-decreasing, conditioned on a stronger uniqueness result for (|2.4p which we 
conjecture to be true. 

From now on, given a probability measure v in flzr {^ex resp.) we denote by F^, 
the distribution of the process ^" (77" resp.) starting from the initial distribution v 
and we denote by Ei^ the expectation with respect to P,y. For the invariant measures 
^p, we wiU write P'' = P^^ and E'' = E^^. 

2.6. Entropy estimates. In this Section we discuss the relevance of hypothesis 
(H). The main point is that the entropy H{i'"{t)\i'p) is decreasing in time, and 
therefore it can be used as a Lyapunov function for the evolution of ■ The results 
in this section are standard; they were introduced in [5] and we include them here 
for the sake of completeness. We will follow the exposition of [9]. 

First, we recall a variational formula for the relative entropy H{p\i'): 



H{p\u) = sup I j fdp, - log j e^di/j, 



/ 

where the supremum is over all functions / which are integrable with repect to p. 
A very useful way to estimate the integral of a function / with respect to p in terms 
of the relative entropy H{p\i') is obtained taking 7/, 7 > as a test function in 
the formula above: 

/d^<£(/M + ilog/e^/di.. (2.8) 
7 1 J 

This inequality is known as the entropy inequality. Of course, this inequality is 
not useful unless we have a good way to estimate relative entropies. For t > 
define /" = dv'^{t)/dvp. The density /" satisfies the Kolmogorov equation 

Notice that in our case the measure Vp is reversible and therefore L*^^ = Lzr, 
although this point is not crucial. Define Hn{t) = H [v"' {t)\v p) . We see that 

d 



dt 



,f„(o^|//r.og/r..,^/(i+.og/f),.»L:../f^.. 

= j n' n' L, ,{1 + log f^jdv, = n' j i;L,Aotfrd„, 
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The operator Lzr, being the generator of a particle system, is of the form 
X^ie/ ~ fid] for some set of indices i € /, some non-negative rates a 
and some transformations ^ i— > ^* of the space VL^r- Using the elementary inequality 
a(log6 — loga) < 2y/a{^/h— y/a), on each one of the terms composing /"L^,- log/", 
we obtain 

Define V{f) = -J VfL.r^diyp. Then, 

|iJ„(0 < -2n"I?(/r). 
Integrating the previous inequality between and t, we obtain that 

H,,{t) + 2n" / V{f^)ds < ff„(0) < Kn'', 
Jo 

where the last inequality is due to hypothesis (H). In other words, entropy is 
decreasing in time, and moreover it also controls the growth of the so-called Dirichlet 
form T>{fl^). Define /" = t^^ f^ds. By convexity of T>, we conclude that 

V{fn < (2.9) 

We can think about H {v^ {t)\v p) / nf^ as a measure of the macroscopic entropy of 
the system. What is remarkable is that the bound (|2.9p on the Dirichlet form is 
enough to control the space-time fluctuations of the density of particles, as we will 
see in the following sections. This observation was introduced in [S] and it is at the 
heart of the proof of hydrodynamic limits for particle systems. 

3. The exclusion process with long jumps 

In this section we prove Theorem 12.31 Let us define the empirical process tt" by 
7rr(dx)=n-'^ ^ 7^l\z)6,/Adx), 

which turns out to be a Markov process in I?([0, oo), A^_|.(R'')). Notice that the 
weak topology in + is metrizable. Theorem l2.3l is an immediate consequence 
of the following result: 

Theorem 3.1. Under the hypothesis of Theorem \2.!A the process {7r";i e [0,T]} 
converges in distribution to the deterministic trajectory u{x,t)dx, where u{x,t) is 
the solution of i2.1\) . 

A standard continuation argument shows that this theorem also holds for the 
process tt" defined in r>([0, oo), M+{W-)). We restrict ourselves to a bounded inter- 
val to simplify some of the arguments. The proof of this theorem follows the usual 
approach to convergence in distribution of stochastic processes. First we prove 
tightness of the distributions of {7r"}„. Then we prove uniqueness of the possible 
limiting points. Since a relatively compact sequence on a metrizable space with 
only one accumulation point is necessarily convergent, Theorem 13.11 follows from 
these two affirmations. 



12 



MILTON JARA 



s2.1 

3.1. Tightness. We recall that a possible choice for a metric that generates the 
weak topology in j\4+(M'^) is 



c?(7ri,'7r2) = ^^minj J Gid(7ri - 7r2), l|, 

i>l 



where {Gi}i is a numerable collection of suitable non-negative functions in C^(M'^), 
the set of infinitely differentiable functions with compact support. Denote 7r"(G) = 
/ GdTr". It is not difficult to see that the sequence {7r"}„ is tight if and only if the 
sequence {7r"(G')} is tight for any non-negative function G G C^{R'^). Notice that 
the projections 7r"(G) are real- valued, and therefore easier to handle than the full 
process tt". By Dynkin's formula, 

7rnG)=<(G)+/ <(£„G)ds + Mr(G), (3.1) 
Jo 

where the operator is defined by 

CnG{x/n) = ^ nXz) {G{x/n + z/n) - G{x/n)^ (3.2) |ec4.0 

and M.^{G) is a martingale. The martingale A^"(G) has mean zero and quadratic 
variation 

{M^G)) = ^ /* -4 E "X-^ - - <(y))'(G(-^/n) - G{y/n)fds. 

We will make repeated use of the identity p{z) = h{z) = n~^'^^°'^h{z/n). Notice 
that 



(A4r(G))<^ J2 ^h{z/n-y/n){G{z/n)-G{y/n)) 



2 



This last sum is nothing but a Riemann sum for the energy £{G, G). A Taylor 
expansion of G shows that the integral defining £{G,G) is absolutely convergent, 
as well as the Riemann sum above. In particular, we have (A1"(G)) < C{G)t/n'^. 
Therefore, the martingale A1"(G) converges to in £^(P^n). In the same way, 

CnG{x/n) = \Y^ h{z/n)(G{z/n - x/n) - G{x/n)^ 

which is a Riemann sum for CG{x/n). Since G G C^iW^) and due to the symmetry 
of h, it is not diSicult to show that 

lim sup I £„G(^/n) -£G(^/n) I =0, (3.3) | ec4.1 

lim n-'^ V \CnG{z/n) - CG{z/n)\ = 0. (3.4) |ec4.2 

The simplest way to prove tightness of the sequence {7r"(G)}„ is to use Aldous' 
criterion, which now we explain. 

Proposition 3.2 (Aldous' criterion). Let (£,d) he a separable, complete metric 
space. Let {P"} be a sequence of probability measures in X'([0, oo), £). The sequence 
{P"}„ is tight if: 
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i) for any e > there exists a compact set K C £ such that 

supP"(7ro i K) < e, 



n 



ii) for any e > and any T > 0, 

limlimsup sup F"(d(7r(T-_,.^-)/^j-, tTt-) > e) = 0, 



"5— *0 Tt— ►oo t£Tt 



7<5 

where 7t is the set of stopping times bounded by T . 

In our case, condition i) is automatically satisfied due to the fact that {/^"}n is 
associated to the bounded profile uq. Condition ii) follows from equation (|3.ip . In 
fact, since the number of particles per site is bounded by 1, 



It is easy to see that |£„G(a;)| < C(G)/(1 + ||a;||'^+") for a constant C{G) de- 
pending only on ||G||oo, ||G"||oo and the support of G (here G" denotes the Hessian 
of G). Plugging this bound into the inequality (|3.5p . condition ii) follows for the 
integral part. By the optional stopping theorem and Tchebyshev's inequality we 
have 



P^.(|A^f,+^)^j,(G) - X;'(G)| >e)< e-^E^,. [(A^^,^^)^^(G) - X;'(G))'] 

^ 6-2E^„ [{Ml+,)^r{G)) - (M"(G))] 
^ g(G)7 



which goes to as n ^ oo, uniformly in r and 7 < i5. In the last line above, the 
constant C(G) depends only on ||VG||oo and the support of G. Here and in the 
sequel we denote by G a generic constant which may change from line to line, but 
depends only on the parameters indicated (for example, in the lines above G(G) 
depends only on G). Therefore, the three terms on the right-hand side of (|3.ip are 
tight, from where tightness for {7r"(G)}„ (and in consequence for {7r"}„) follows. 

3.2. Uniqueness of limit points. Once we have proved tightness for {tt"}, we 
know that this sequence has accumulation points with respect to the topology of 
convergence in distribution. Let tt. be one of these points. Denote by n' a subse- 
quence for which tt" converges to tt. The idea is to to take the limit through the 
subsequence n' in p.ip . By definition, tt" (G) converges to iTt{G). By the assump- 
tions on initial distributions, ttq (G) converges to J G{x)uo{x)dx. The martingale 
term A4"{G) converges to in £^(P^^) and in particular it converges to also 
in distribution. However, £„G is not a function in Cc(K''), due to the non-local 
character of the operator £„, so some justification is needed before taking the limit 
through n' of 7r"(£„G). Observe that the number of particles per site is bounded 
by 1. Using p.3p . we can substitute CnG by £G in (|3.ip by introducing an error 
term that vanishes as n —^ 00. Notice that 




(3.5) 



ec5 




\z\>Mn 
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uniformly in n. In particular, we can approximate CG by functions of compact 
support to obtain that 



lim / <'(/:G)ds= / TTs{Cg)ds, 
Jo Jo 



in distribution. Therefore, the measure tt. satisfies 

Jo 

for any function G G C^{R'^). Let G* : [0, T] C^{W^) be a differentiable trajec- 
tory. Repeating the arguments above, it is not difficult to prove that tt. satisfies 



^t(Gt) = 7ro(Go) + / TTs{{dt+C)G)ds. (3.6) 
Jo 

Since the number of particles per site is bounded by 1, it is not hard to see that 
the limiting measure ttj is, for any t > 0, absolutely continuous with respect to 
Lebesgue measure. Moreover, its density is bounded between and 1. Let us write 
TTtidx) — u{x,t)dx. In terms of the random (at this point) density u{x,t), equation 
p.6p states that 

u[x,t)Gt{x)dx — j uo{x)Go{x)dx + j j u{x, s)(^dt + £)G{x)dxds 

for any smooth trajectory Gt, which is exactly the weak formulation of the hydro- 
dynamic equation (|2.ip . In other words, we have proved that tt. (dec) = u{x,-)dx 
is concentrated on weak solutions of (|2.ip . But this equation has only one weak 
solution starting from uq (see Section [8.1|) . This uniqueness result finishes the proof 
of Theorem 13.11 



4. The zero-range process with long jumps 

s3 I 

As we did for the exclusion process, we will consider the empirical measure 

<idx) - ^ E ^Tiz)Sz/nidx) 

and we will prove that 

t4 I Theorem 4.1. Under the hypothesis of Theorem \ 2.5l the sequence {7r"}n is rel- 
atively compact with respect to the topology of convergence in distribution in the 
Skorohod space X>([0,r],7W + (M'^)). All the limit points are concentrated on finite 
entropy solutions of the hydrodynamic equation (2^. When \2.4-\j has a unique 
solution u{x,t), tt" converges in probability to the deterministic path u(x, ■)dx. 

We will explain what we mean by a finite entropy solution in Section [71 where 
we state a more precise version of this Theorem. The proof of this Theorem follows 
the same strategy followed in order to prove Theorem l3.1i but it is technically more 
involved. For a proof in the case of a finite range, mean zero transition rate p{-), we 
refer to Chapter 5 of [9]. For the reader's convenience, we follow closely the proof 
in [9], modifying to our setting. 
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s3.1 I 

4.1. Some elementary estimates. Before we enter into the proof of Theorem l3 . 1 [ 
in this section we collect some elementary estimates that will be used repeatedly. 
The results are well known, and we collect them here for the reader's convenience. 
Let G G C~(R'^) be given. Define ^0(2:) = 1/(1 + Our first estimate gives 

the behavior of £„G and CG at infinity. 

auxl Lemma 4.2. For any function G G C^(M'^), there exists a constant C{G) such 
that for any 2 G and any x G M'^, 

\CnG{z/n)\<G{G)Fo{z/n), 

\CG{x)\ < G{G)Fo{x). 

Proof. Observe that, since h{-) is strictly positive, there exists a constant eo such 
that eo||x||-(''+") < h{x) < for any x ^ 0. Therefore, it is enough 

to consider the case h{x) = For this last case, the result of the lemma 

follows easily from the compactness of the support of G and a Taylor expansion of 
second order. □ 

The second estimate concerns the behavior of the moment generating function. 

aux2 Lemma 4.3. The function log M p{-) is strictly convex and increasing. In particu- 
lar, if Mp{9o) is finite, 

logMpie) < 



Mp{Oo) 



70 

for any < 0o- 

Proof. Let us denote by qp.e the distribution in No with density e^" /Mp{9) with 
respect to <?p. Let us denote by Ee the expectation with respect to qp,e. It is enough 
to observe that for ^{9) = log Mp{0), ^''(0) = Ee[n] and ^"(e) = Vare[n]. □ 

For a given function G : R"* ^ M, let us define 

aux3 Lemma 4.4. Let G : M'* M. be a bounded function. There exists a constant 
G > such that for any t > we have 

E.4|7rr(G)|] <G(||G|U + ||G||i,„). 

Proof. Without loss of generality, we can assume that G is non-negative. Remember 
that E„r.[nl'{G)] = J TT"-{G)f^dvp. By the entropy inequahty. 



E..K(G)] < ^ + -Liogl exp{7 ^ ^{z)G{z/n)]dvp 



K 1 

7 7n" 



<- + —dl^ logAf,(7G(z/n)). 



Take 0q> Q such that Mp{0o) is finite. Then take 7 — 0o/||G||oo- Since the measure 
Vp is of product form, using Lemma [43l we can bound E^™[7r"(G)] by 

K\\G\U \\G\U Mpi0o)G{z/n) 

Oo 0on<^ f^, IIGIloo ■ 

Taking G = 0'^ maxjif, Mp{0Q)} we obtain the desired bound. □ 
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4.2. Tightness. Remember that in order to prove tiglitness for {7r"}„, it is eiiougli 
to prove tightness for {7r"(G)}„ for any non-negative G £ C^(M'*). The projection 
of the empirical measure can be written as 

7rr(G) = <(G) + r ^ E 9{Ci^))^nG{z/n)ds + MnG), (4.1) 

where A^"(G) is a martingale of quadratic variation 

{MnC)) = T-L ^ giCiy)) E h(z/n-y/n){G{z/n)-G{y/n)yds. 

Let us define 

QnG{x/n) ^^^/"^ ~ x/n){G{z/n) ~ G{x/n)f. 

Notice that QnG{x/n) is a Riemann sum for the integral / h{y — x){G(y) — 
G{x))^dy evaluated at x/n. We can rewrite the quadratic variation of A^"(G) as 



(MnG)) 



Ef(eW)QnG(z/n)ds. 



n 



In this expression, as for the exclusion process with long jumps, there is an extra 
1/n^ term that will make the quadratic variation of Ai"{G) converge to 0. Now 
the difference is that we need to bound the expectation of g{^l''{x)), which is not 
longer a bounded random variable. We can write (A^"(G)) in terms of 

(MnG)) = ^ 9{i{z))Qr.G{z/n)r;dv,. 

Remember the bound g{n) < nn. Notice that the entropy of /" is also bounded 
by Kn'^. Therefore, the proof of Lemma [4.41 also applies here, so we obtain the 
bound 

Cf 

m,^{M\\G)) < ^(||Q„G||oo + ||Q„G||i,„). 

Notice that both norms ||Q„G||oo + ||QnG||i.„ are uniformly bounded in n. We 
conclude that A^"(G) converges to in £^(P,yii) as n goes to oo, like in the case of 
the exclusion process. Therefore, {A1"(G)} is tight. 

With the same notation of Section 13.11 and using the entropy estimate, for any 
non-negative, bounded function F we have 



(t+7)AT 



J2 giCi^m^Mds 



K 
< — 

70 

f.(r+7)AT 



exp {^oJJ ^ E 9iCi^))Fi^/n)ds} 



HYDRODYNAMIC LIMIT OF PARTICLE SYSTEMS WITH LONG JUMPS 



17 



By Jensen's inequality, this last line is bounded by 

(t+7)AT 

< - + ^\ogEp\- [ exp|7o7 ^ 9iCiz)))F{z/n)}ds 
70 70'^'' L7 J I ^ J 

logE^ji rexp{7o7 ^ 



~ 70 70'^'' 



where we have used the fact that g{n) < nn in the last line. Since the measure Up 
is invariant under the evolution of the expectation is bounded by 



K 

70 70"" '^7 



i^logl- J] Mp(7o7/«F(z/n))| = 

{log(T/7)+ E log Mp{^o^KF{z I n))y 



70 7o?^" 

For 7 small enough (how small depends only on ||i^||oo)j by Lemma 14.31 the expres- 
sion above in bounded by 

- + -^A log(r/7) + E C{p){e^U^nF{z/n)\, 

where the last inequality is true for 7 small enough and Ci^p) si simply Mp{0q)/0o. 
We conclude that 

(r+7)AT 

sup limsupP„,.f / ^ V 9{C{z))F{z/n)ds > e\ < — + C{p)5n [ F(x)dx 
7<5 n^oo V 7 71" ^ 70 J 



7 



ze 



for 6 small enough, if F is such that the Riemann sum above converges as n — s- 00. 
For G e C^{M.'^) and using Lemma [4.21 tightness follows for the integral term in 
(|4.ip by Aldous' criterion. This finishes the proof of tightness for the sequence 
{^"(G)}„. 

4.3. Identification of limit points. Let G :W^x [0, 00) ^ R be a function of class 
C^'^ and bounded support. Write Gt{-) = G{-,t). Following the same computations 
leading to tightness made on the previous section, we see that 

TrHGt) =^S'(Go)+ / 7r^{dtGs)ds+ 
Jo 

+ [ ^J29iCi^))^nGsiz/n)ds + MnG), (4.2) 

^0 " za^ 

where Mt{G) is a martingale such that E^n [7U"(G)^] — > as n ~> 00. We already 
know that the sequence {7r"}„ is tight, so now we want to characterize its limit 
points as energy solutions of the hydrodynamic equation (12.41) . The main difference 
between p. II) and (|4.2p is the presence of the function g{^"{x)). This makes the 
second integral in (|4.2p not to be a function of the empirical measure tt". Our 
objective, then, is to write the integral 

^ E 9iC{^))^nGs{z/n)ds 
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as a function of the empirical measure plus a rest that vanishes in distribution as 
n ~f oo. This is the content of the following theorem, known in the literature as 
the Replacement Lemma: 

Theorem 4.5 (Replacement Lemma). Under the hypothesis of Theorem \ 2.5[ for 
any M > and any t > 0, 

limlimsupE,. / — ^ \ ^ g{C{x + y)) 



\x\<AIn |y|;^€n 



\y\<en 



ds 



Here and below, by abuse of notation we write en in place of its integer part 
[en] = sup{n e Z; n < en}. Notice that with this convention, (2en + 1)'^ is just the 
cardinaHty of the set {\y\ < en}. We postpone the proof of this Theorem to the 
next section. The Replacement Lemma is, with no doubt, the heart of the proof of 
Theorem l4.1l Let tt. be a limit point of {7r"}n, and denote by n' a subsequence for 
which tt" TT. in distribution. Notice that 

(^2en+ir ^ Cix + z) = c^n-{{2er'li\y^x/n\<e)), 

\z\<en 

where c„ is a normalizing constant that goes to 1 as n ^ oo. In particular, this 
expression is a function of the empirical measure. Define, for any measure tt G 
M + {R'^), the function J.vr by 



IM^) = (2e)-^ / nidy). 

J \y—x\ <e 

Since the function £„G's(x) is uniformly continuous and bounded by C{G)Fo{x), 
thanks to the Replacement Lemma we can write 



f \Y.9{Cs{^))'^nG{x/n,s)ds 



* 1 V- 



= / —2^Hlen:{x/n))C,,G{x/n,s)ds + R'^'%G), 

where R"''^{G) is an error term that vanishes in C^(¥,^n) when n — > oo and then 
e ^ 0. We denote by R^''^{G) any term with this property. 
Using (|3.3p and p.4p . we can write 

— Yl H^eKix/n))CnG{x/n,s)ds = 

= / / 4>{leTT'^{x))CG{x,s)dxds + R'/'iG). 
Jo Jw^ 

Notice that we can not say that /eTt" (x) —> I^tts^x), since the indicator function 
1(1 ■ < e) is not a continuous function. However, by Portmanteau's lemma, this 
is true whenever 7rs(9{t/; |t/ — ccj < e}) = 0, where we denote by dA the boundary 
oi A C M."^. Since tt is a Radon measure, this is the case in a set of full measure 
in K'^. In particular, I^tt" (x) IeTTs{x) for a; in a set of full measure in R"*. 
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Since \(j)[len^{x))£G{x,s)\ is bounded by KleTT^{x)\CG{x, s)\, by the dominated 
convergence theorem we conclude that 

lim / / (jyll^TT"^' {x))CG{x,s)dxds = / / (j)(le'Ks{x))CG{x,s)dxds. 

Therefore, taking the limit in (|4.2p through the subsequence n' , we obtain that 

G{x,t)TTt{dx) ^ I uo{x)G{x,0)dx + / / dtG{x, s)TTs{dx)ds 



(j>{l,TTs{x))CG{x, s)dxds + RUG), 

where Rl{G) is a rest that vanishes in expectation when e — > 0. At this point, we 
only need to find the limit as e ^ of 

ft f 

0(/e7rs(a;))G'(x, s)dxds. 

It is clear that is just an approximation of the identity, so it is reasonable to 
expect that I^-Ks — > tt^ as e — > in some sense. We need some regularity for tTs in 
order to pass to the limit inside the function 0. This is the case if, for example, tt^ 
is absolutely continuous with respect to Lebesgue measure. The following lemma 
says that this is indeed the case. 

1 11 I Lemma 4.6. The process n. is concentrated in trajectories of the form u{x, ■)dx, 
where u : [0, T] x — > M is non-negative and locally integrable. Moreover, 

•-T I , , 

H{u{x, t))dxdt < +00. (4.3) ec8.2 



We will asume this lemma and we will prove it later. Since TTs{dx) = u{x,s)dx 
and u(x,s) is locally integrable, we see that I^t:s(x) converges to u(x,s) a.s. with 
respect to Lebesgue measure, and also in £l^^{R'^). The bound g{n) < Kn implies 
that (j){p) < Kp. This linear bound plus the fact that (f){p) is locally Lipschitz (since 
it is smooth) allow to conclude that (j){IgTTs{x))dx converges weakly to (j){u{x, s))dx 
as e — + 0. By Fubini-ToncUi's theorem, we conclude that 



lirn j j (t)(^I^TTs{x))CG{x, s)dxds = j j 



(f>{u{x, s))dxds. 



Therefore, u{x, t) satisfies conditions i) and iii) of the definition of energy solu- 
tions of (|2.4p . In order to finish the proof of Theorem 14.11 it is only left to prove 
that u{x,t) satisfies the energy estimate J^^ £{(f){ut), (l){ut))dt < +oo. This is the 
content of the following theorem: 

t5.1 Theorem 4.7 (Energy estimate). Under the hypothesis of Theorem \2.5[ the limit 
points of {7r"}„ are concentrated in measures of the form u{x, ■)dx, where «*(•) = 
u(-,t) satisfies 

rT 

£{4>{ut), 4>{ut))dt < +O0. (4.4) I ec8.1 

Accepting the validity of Theorems 14. 5i 14.71 and Lemma 14. 6i we conclude that 
the limit point tt. is concentrated on enegy solutions of the hydrodynamic equation 
(|2.4|) . We will dedicate the next sections to the proof of Theorems 14.51 14771 We end 
this section proving Lemma [4.61 
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Proof of Lemma The reader may recognize in the following lines an application 
of Varadhan's lemma. Let us define h{9) — \ogMp{6). For any sequence of random 
variables ai , . . . , we have the estimate 

logi?[exp{max{ai}}] — log i?[max{exp{ai}}] < logi?[y^ exp{ai}] 

i 

= log^i?[exp{aj}] < log {fcmaxiJ[exp{ai}]} (4.5) 

i 

< log k + max log £^[exp{ai}] . 

i 

Let , . . . , G'"' be a sequence of functions in Cc(M'' x [0,T]). We will apply 
the estimate above to the sequence /g^ ■K^{G\)dt, ■ ■ ■ , Jq iT"{G'j;)dt. For a function 
G :R'' ^R, let us define J„(G) = 'n'~'^J2x h{G{x/n)). By the entropy inequality, 



max 

i— 



< 



< X + — rlogEP 

n"- 



cxp { ^ max^ ^ 1^ (7r"(GD - Jn{Gl))dt} 



< K ^+ max — rlogE' 

n" i=i,...,k n" 

log A: 1 , 1 

< A H — h max —7 log 



i=i,...,fc n" 



exp ■ 

T 

E'' 



{y I {^''{G\)-UG\))dt] 







exp{n''(7r"(Gj)-J„(Gj))} 



dt. 



ec9 



Taking the limit along the subsequence n' , the left-hand side of this inequality 
converges to ii^[maxi Jo {'^tiGD — J{Gl)}dt], and the right-hand side converges 
to K. Now k is arbitrary, so we have proved that 

e\sup f {Trt{Gt)-J(Gt)}dt\ <KT, (4.6) 

where the supremum is over G G Cc,t —'■ Cc(M'' x [0, T]). The supremum inside the 
integral can be computed. Denote by TC the Legendre transform of h: Ti.{a) = 
supg {a9 — h{9)} . By Lemma [4.31 h{-) is strictly convex and smooth. Therefore, 
7i(-) is also strictly convex, and grows at least linearly. Let us define 

J{tt)= sup / {7rt{Gt)-J{Gt)}dt. 

Gec,,T Jo 

Suppose that ttj (dx) is absolutely continuous with respect to Lebesgue measure. 
Let u(x,t) be its density. Then, 

Jin)^ sup / I {u[x,t)G{x) -h{G{x))]dxdt 

G<^Ca,T Jo J 

< / snp \^u{x,t)g — h{g)}dxdt ^ / / Tl{u{x,t))dxdt. 

Jo J geR Jo J 

By Lemma |4.3[ h{-) is strictly convex and smooth. Therefore, in the definition of 
7i(a) the minimizer 9{a) such that 7i(a) = a9{a) — h{0{a)) is unique and smooth as 
a function of a. In fact 9 = 9{a) satisfies a = h'{9). After some computations, we 
can see that the function Ti.{a) is equal to the entropy function defined in p.Sp . We 
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conclude that the supremum in the definition of i7(7r) is attained at the continuous 
function 6{u{x,t)), and therefore J^{tt) = j Ti.{u{x,t))dx. Now assume that ttj is 
not absolutely continuous with respect to Lebesgue measure. We will prove that 
1/(71) = +00. Let F be a bounded, closed set such that ■Kt{F)dt = a > and 
Jpdxdt = 0. Fix /3 > 0. Let F, = {{y,t) e R'^ x [0 , T]; d{{y , t) , F) < e} be the 
ball of radius e and center F. Let : R'^ x [0, T] [0, (3] be a continuous function 
such that Cix, <) = /3 if {x, t) e F,/2 and G'{x, i) = if (x, t) F,. Then, 

nt{Gl)dt- [ J{Gl)dt>aP-h{f3)[ dxdt. 
Jo Jf, 

Sending e — s- 0, we conclude that i/(7r) > a/3. Since P is arbitrary, we deduce that 
J{Tr) = +00. By (|4.6|) . we know that E[J'{it)] < K. In particular, J{t^^ is finite 
a.s., from where tt. is concentrated in trajectories of the form u(x, ■)dx. Since h{-) is 
strictly convex and smooth, 7i(-) is also strictly convex, and grows at least linearly. 
We conclude that u{x,t) is locally integrable. □ 



5. The Replacement Lemma 

In this section we prove Theorem 14. 5 1 Following [5J, the proof is divided into two 
pieces: the so-called one-block and two-blocks estimates. The one-block estimate 
does not pose a real challenge for the model we are considering. After proving 
the one-block estimate, the two-blocks estimate follows from the moving particle 
lemma, which roughly states that the cost of moving a particle from one site x to 
another site y can be estimated by the Dirichlet form of the process restricted to 
a box containing both points. In the diffusive case, the moving particle lemma is 
an elementary application of Cauchy-Schwartz inequality, but in our superdiffusive 
setting, a more sophisticated proof is needed. For the sake of completeness, we give 
a proof of both one-block and two-blocks estimates. 

To simplify the notation, for a; e Z'' and / > we define 
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y), 



\y\<i- 



In terms of , the Replacement Lemma can be written as 



lim lim sup Ej^„ 

^ ^0 n — voo 



= 0. 



(5.1) 



\x\<Mn 

The expectation in (|5.ip can be written as 



\x\<Mn 



It will be convenient to introduce an intermediate scale / into (|5.ip . The limit 
in (|5.ip is a simple consequence of the following two lemmas (see [9], Sect. 5.3): 
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I lb I Lemma 5.1 (1-block). Under the hypothesis of Theorem 



lim limsup / / — 7 

n^oo J J 

\x\<Mn 



2b Lemma 5.2 (2-blocks). Under the hypothesis of Theorem \2.5\ 

lim lim sup lim sup sup [ + z) - ^'■{x + y)\f"{£_)di^p ^ 0, (5.2) 



where the supremum is over y^z ^1/^ such that |?/|, \z\ < en, \z ~ y\ > 21 + 1. 

Observe that "pasting" boxes of size n in a convenient way, we can assume 
without loss of generality that AI = 1. Therefore, from now on we assume M ^ 1. 



5.1. The one-block estimate. In this section we prove Lemma [5.11 The proof 
is standard and we basically repeat the proof in Section 5.4 of [9J. We include a 
complete proof in order to simplify the exposition of the proof of Lemma 15. 2[ for 
which new arguments are needed. 

Remember that the density /" satisfies the estimates J /"log/"(it'p < Kn'^, 
2?(/r) < Kn'^-°'/2t and fj'vp di V- According to Section 5 of [9 , this is all we 
need to know about In fact, the coupling estimate /"j^p di i^p' is only needed to 
prove the two-blocks estimate when g{-) satisfies (B). First we introduce a cut-off 
of large densities. As in step 1, Section 5.4 of [9 when g{-) satisfies (B) or Lemma 
5.4.2 of P when g{-) satisfies (FEM), the following limit holds: 

lim limsup limsup I \ V^{S,)1{S} {x) > a}f1'dvp ^ 0. 

n — *oo J n II 

\x\<n 



Let us define Vf:'°'{^) = T4'(^)l{^'(x) < a}. In order to prove Lemma [57T1 we just 
need to prove that for any a > 0, 

1 . 



lim sup lim sup 



I a: I < n 



(5.3) 



From now on we simply write ft=f, keeping in mind the properties of /" stated 
above. Let denote the translation by x: TxS,{z) = £^{x -t- z). Performing the 
change of variables ^ = TxS, and by translation invariance of the measure Up^ we 
can pass the sum in x from to /. In this way we see that the integral in ()5.3p is 
equal to 



1 



{2n+ ly 



(5.4) 



where we write TxfiCi = f{TxCj- The prefactor {{2n + l)/nY is bounded and it can 
be ignored. Notice that the weighted sum in (|5.4p is also a density. We denote this 
density by /. Define A/ = {— /, . . . ,1}'^ and let Ti be the ti-algebra generated by 
{^(x); X e A;}. For simplicity, we also denote by A; the cardinality of {— /, . . . , 
which is equal to {2l + \Y. Notice that A„ = {x G Z''; |a;| < n}. Define /; = E[f\Ti], 
that is, // is the conditional expectation of / with respect to the configuration ^ 
restricted to the box A;. With this notation, and due to the product structure of 



the integral in (|5.4p is equal to 



vtmrn^p, 



(5.5) 



ecl2.1 



ecl3 



ecl4 



ecl5 
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where now the integral is over the set fli — Ng ' . At this point we have reduced 
the original problem into a finite-dimensional, static problem. In fact, due to the 
indicator function, the integrand is difi'erent from zero only on the finite set G 
f],;C'(0) <a}. 

The following step is to estimate the Dirichlet form of //. For any e Z*^, let 
us define 

2 



dvp. 



Notice that V'y{f) = X>^'"^(/), due to the symmetry of p(-). The form V^vif) 
is convex as a function of /. Since / is an average over translations of /, we have 



For a given density /, let us define 

The form !?;(•) corresponds to the Dirichlet form associated to the generator L^,^ 
of the process restricted to the box A;. We have 

x,yeAi 
2GA„ 

By the convexity of T>i{-) is also convex and Vi{fi) < !?/(/). We conclude 

that 'Di{fi) < c{l)/n°', where c{l) is a constant that only depends on /, K and t. 
Therefore, the integral in (|5.5[) can be estimated by 

sup / V^iOfiOdiyp. (5.6) 

■Di(/)<c(i)/n""' 

Due to the indicator function in the definition of Vq'"", we can restrict the supre- 
mum to densities supported on {^'(0) < a}. This set of densities is compact, and 
both the function 'Di{f) and the integral in (|5.6p are continuous. Therefore, the 
limiit as n — > oo of (j5.6[) is bounded above by 



sup / K)'"(0/(Orf^P- (5.7) 
v,u)=oJ 

If 'Di{f) — 0, then the density / is constant on each ergodic component of the 
state space ili with respect to the dynamics generated by L^,^. Since the number of 
particles is the only conserved quantity with respect to this dynamics, these ergodic 
components are exactly the sets S f2;;f'(0) — k}, k > 0. Denote by Vk,i{') the 
measure Vp{-\S}{Q) = k/Ai). Then, the limit as Z ^ oo of the supremum in l|5.7p is 
equal to 



lim sup / Vo{^)di^k,i 

fc<aAi J 



l^oo k<aAi 

This last line is exactly the same limit appearing in pg. 89 of |9j, and it is equal 
to due to the equivalence of ensembles. This ends the proof of Lemma [5?T] 
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5.2. The 2-blocks estimate. The proof of Lcmma l5.2l is similar to the proof of the 
one-block estimate explained in the previous section. When g(-) satisfies (FEM) 
(see section 5.5 of [5]), using the entropy inequality we can prove that 

lim lim sup lim sup sup / — |^'(a; + z) — ^'(x + t/)| x 

\x\<AIn 

X l{max{^'{x + y),^'{x + z)}>a}fl^{Odiyp = 0. (5.8) 

When g(-) satisfies (B), we need to use hypothesis (C) to introduce this cut-off. 
The integrand in (|5.8p is increasing and f^i^p ^ Vpi. Therefore, we can estimate 
the integral replacing /"(^)rfi^p by dvp'. By Tchebyshev's inequality, the integral 
can be estimated by 4 / ^{x)^di'p' /a, which proves that (|5.8p also holds when g(-) 
satisfies (B). Define ^^,^(0 = < a}l{?'(y) < «}■ As we did in the proof 

of the one-block estimate, passing the sum to the density /, the integral in (|5.2p is 
equal to 

[-^) J \e{z)~e{y)\BiMm>p{o- (5.9) 

Notice that the function l^^z) — £,\y)\ depends on the configuration of particles 
inside the two disjoint blocks TyAi and TzAi. The name "two-blocks estimate" 
comes from this observation. Let us denote by A^* the union of these two blocks, 
and let us write ^i(a;) = ^{x + y), £,2{x) — £,{x + z). We represent A;* as A; x A;, 
dropping y and z from the notation. The integral in (|5.9p is equal to 



A; xA, 



l^i(o) -^^(O)|sr(0ii*(ei,6)^p(d6 X d6), (5.10) 



where /;* = E[f\T\*], Tiy* is the cr-algebra generated by {^(a;);a; € A^*} and 
B^i) = 1{?1(0) < a}lU^(0) < a}. 

The integral in (|5.10p depends on y, z only through ff. The estimation of this 
integral is now a finite-dimensional problem. This integral is similar to the one in 
(|5.5p . so we will estimate a suitable version of the Dirichlet form of /;*. Let us 
define 



\x\,\x'\<l 

p(z - y) 



Following the computations made when we estimated !';(/;) in the previous 
section, we see that 

vioifn < (5.11) 

The term I?*''^(/;*) is the one that connects the behavior of /;* between the two 
blocks, so it is the most relevant in the definition above. Notice that the dynamics 
associated to !';*(•) corresponds to a system on which particles perform a zero- 
range process with long jumps restricted to each one of the two boxes, and on 
which particles can jump from site y to site z (site z to site y resp.) with rate 
g{£,{y)) igi£.iz)) resp.) The rate at which particles jump between y and z does not 
depend on the distance between y and z. Now we are ready to state what we mean 
by the moving particle lemma: 
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Lemma 5.3 (Moving particle lemma). Under the hypothesis of Theorem \2.5\ 

limlimsupsup— — zVy-^'ifr) = 0, (5.12) 

e^O n^oo y,z P(Z - y) 

where the supremum is over y, z € A^n such that \z — y\ > 21 + 1. 

If we understand by p{z — ?/)~^2?^'^(/*) as the cost of moving a particle from y 
to z, what this lemma is saying is that the cost of moving a particle at a macro- 
scopically small distance vanishes when the distance goes to 0. We will prove this 
lemma in the next section. Assuming this lemma, it is not difficult to finish the 
proof of Lemma [Ol By (|5.1ip and (|5.12p . X',*(/;*) goes to as 7i — > cxj and then 
e ^ 0. Therefore, by the same compactness argument used in the proof of the 
one-block estimate, we can bound the limit of the integral in (|5.10p by 

hm sup / |^kO)-e^(0)|i3r(eX,(da xd6), 

fc<2a/A, J A, xA, 

where i'^ ^ is the measure in A/ x A/, conditioned to have exactly k particles. 
This limit is equal to due to the equivalence of ensembles, which ends the proof 
of Lemma 15.21 



6. The moving particle lemma 



In this section we prove Lemma l5.3l To avoid heavy notation, we start assuming 
d — I. Later we explain how to generalize the proof to arbitrary dimensions. 
Assume, without loss of generality, that y < z. Let {y = y'^,y^, . . . ,y"^ = 2;} be a 
path from y to z. Notice that 



e=^(((..v)''-)...)' 



In other words, moving a particle from y to 2: is the same that moving a particle 
from y to y^ , then from y^ to y^, etc. Using the formula above, write 

m 

1=1 

Using the inequality (ai + • • • + 0^)^ < m(a^ + • • - + 0^,^) and performing a change 
of variables, we conclude that 

— -^<r(iy—. ^ (6.1) 

p(z-y) fr'^piy'-y' j 

for any density /. Now the idea is to choose the path between y and z in an 
adequate way. The simplest choice is = y + j. In that case, 

_ |'^' P^+"^-^+'(/) 

p{z-y)-\' p{i) ■ 
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Applying this estimate to /;*, we see that 

i=l,...,\z-y\ 

- Pil)n Pil) ■ 

The last but one inequality follows from the fact that each of the terms of the 
form V^+y '^^y (/) appears at most |z — y| times in the sum. Notice that in the 
diffusive case, a = 2 and this estimate is enough to prove the lemma. The point 
is that in the previous estimate we have only used jumps of lenght one to move 
the particle from y to z. Therefore, we need to use jumps of various lenghts to get 
the right estimate. Fix a positive integer k. Assume for a moment that z — y is 
divisible by k. Defining y^ = y + ik, (|6.ip gives us the estimate 

\x — x^\ — k 

We make use of the following observation: if a^, bk, /3, 7 are non- negative num- 
bers such that ai + • • • + a.,„ < 7 and (3 < h^ak for any fc, then 



/3< 



7 



In our case we take at = T.\x-x'\=k'^'''''' if)^ h = Ae^nk"'\ (3 = Vy^^f^) and 
7 = V{f). Therefore, 

where C{a) is a constant that only depends on a. Here we used the fact that 
Sfc=i k^~°' is of the same order than x^~"dx. Notice that this estimate is still 
not justified, since we assumed that fc is a divisor of z — y for any k between 1 and 
en. When k is not a divisor of z — y, the path from y to z must contain a jump of 
size different from k. The idea is that we can restrict the range of summation of k 
to an interval of the form [aen, ben] , since in that case the sum will still be of order 
(en)^~". Define m — \z — y\. We can assume, without loss of generality, that m is 
divisible by 6. In fact, for Z > 6, what we can do is to take z' at distance at most 5 
from z such that z' — y is divisible by 6, and to consider P^'^ (/) instead of 'Dy'^{f). 
We write m = 6mo. For k = 2rno + Jj i = li • ■ ■ i "Iq, we have m = 2k + 2(mo — j). 
Therefore, we can move a particle from y to z by making two jumps of length k 
and two jumps of lenght mo — j. Notice that k runs from 2mo + 1 to 3mo. 

In this way we have gained control over the rest of the division of m by fc. The 
good property of this decomposition is that for fc ^ fc' we have j 7^ j' and therefore 
we are not repeating jump lenghts for different fc's. Therefore, for any density / we 
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have 



' i=l,...,4 

where y'^ — y, y^ — y + k, y'^ — y + 2k, y^ — z~ {rriQ — j) and = z. Repeating 
the calculations we did to estimate 'Di{fi) in the proof of the 1-block estimate, we 
obtain that 

piz^y) - p{l)n ^ I ^J)-^'^ yj'l 

Since the jumps have different lenghts for different choices of k, we can perform 
a sum from j = \ to j = rriQ to obtain the bound 

P{z -y) ~ p{l)n" - p{l) 

where we have used the estimate mo < en in the last inequality. This ends the 
proof of Lemma 15.31 in dimension d — 1. 

In dimension d > 1, we proceed as follows. First observe that due to the positivity 
and continuity of q{-) on the sphere the Dirichlet form T>{f) is equivalent 

to the one associated to the zero-range process associated to the transition rate 
p{z) = c/||z||''+". In particular, there is a finite and positive constant cq such that 
^o/\\y — x\\ < p{y — x) < ^Q^/Wy — xW^^" To simplify the notation, take d — 2. 
Assume first that (z — y) • 62 — 0, where 62 = (0, 1) is the second element of the 
canonical basis in R"*. Define m = |z — y\. Using the same argument above, we can 
assume that m = Gtoq for some positive integer mo. Define mo = (mo,0) and for 
i, j = 1, . . . ,mo, define k = 2mo + k' = mo — («, j). We can move a particle 

from y to z by making two jumps in direction k and then two jumps in direction k'. 
Now we observe that when i,j run from 1 to mo, the vectors k, k' are all different. 
Repeating the computations done in the d ~ 1 case, 

^' \x\<n{l+e) 

Here C(d, a) is a constant which only depends on a, the dimension d and eq. 
Although we are assuming d = 2, we include the dependence in dimension of this 
estimate for clarity. Performing a sum over i,j e {1, . . . ,mo}, we conclude that 

p(z — y) n" 

When [z — y) ■ ei = 0, we simply exchange the role of the first and second coor- 
dinates and we obtain the same estimate. For y = {yi^y^), z = (zi, Z2) in general 
position, we define y' = (^1,^2). Then, Vy'''{f)/p{z - y) < 2{vy-y' {f)/p{y' - y) + 
'^{f)/p{z — y'), and these two last terms can be estimated as above. 
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7. The energy estimate 

In this section we prove Theorem 14.71 We start introducing some notation. For 
a function F : Z"^ x Z"^ ^ M and a configuration ^ e ^zr, we define 



2 



The proof of Theorem 14.71 consists basically on combining the Replacement 
Lemma 14.51 and the following lemma: 



Lemma 7.1. Under the hypothesis of Theorem \2.5[ for any finite sequence of func- 
tions i^* : [0, T] X Z'^ X Z'^ ^ M, i = 1, . . . , fc such that Fl y(t) = if x,y are big 
enough, we have 



sup 

l<i<fc Jo 



[Ei{i\\F\t))-T^{Q.F^m]dt 



< K 



log k 



(7.1) 



Proof. By the entropy inequality with •f ~ n"^, the expectation in (|4.5p is bounded 

by 



expjn'^ sup r {£^{^^,F\t))~r,,i^-,rit))}dt} 



By estimate (|4.5p . the expression above is bounded by 

T 

f ^ / 

exp 



^_^logfc_^ sup ^logE'' 
i<i<k n<^ 



Therefore, it is enough to prove that 



{n" ^ {^^^(er, F\t)) - r„(er, F^(i))}d4 



-^logE''[exp{n'^ / [s^iS,'^ ,F{t)) -T^{£,'^ ,F{t)]dt] 



< 



(7.2) 



for any trajectory F{t). Since now the expectation is with respect to the process in 
equilibrium, powerful variational methods are available to estimate this expectation. 
By Feynman-Kac formula, the left-hand side of (|7.2p is bounded by n~'^ X{t)dt, 
where X{t) is the largest eigenvalue in £^(fp) of the operator n"_L^'' + n'^V" and V" 
is the multiplication operator given by 

vr/(e) = {m^, F{t)) - r„(^, F{t))}f io. 

By the variational formula of the largest eigenvalue of an operator in >C^(t'p), we 



see that the left-hand side of (|7.2p is bounded by 



1 



sup {{Vr, f^)p - ^{f, -L^'-f)p]dt, 

fdC^iy,) 



^^f^-L'^f)Adt, (7.3) 

where (•, ■)p denotes the inner product in CT'ivp). This integral only involves the 
equilibrium properties of the dynamics. Notice that {f,—L^^f)p = 'D{f^). Let 
us estimate each one of the terms in F{t)), f{£_)^)p separatedly. Define (^ -I- 

Sx){z) = ^(z) -I- Sx{z). Using the change of variables £, ^ £, + Sx, we see that 

/ (g(e(2/))-5(e(a;)))/(e)'d^p = 0(p) / {fii + Syf~f{c + Sx)'}d^p. 
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The same change of variables shows that 

v-^^f ) = Hp) J p{y - + 5y) - f{i + 5..)fdvp. 

Notice the similarity between these two formulas. Using the formula a?' — iF' — 
(a + 6)(a — 6) and the weighted Cauchy-Schwartz inequality 2\ah\ < (3a? + 5^/ /3, we 
see that 

|/(? + 5yf - /(C + = {/(e + 5y) + /(e + 5.)}\f{i + 5y) - /(e + 5^)\ 

<^{/(e+^.)+/(e+4)}' 

+ 7^{f{^^5y)-f{£, + 5^)]\ 

The idea is to choose f3xy in such a way that the second term above cancels with 
n»-dx>^.v(f). This happens choosing fS^^y = \Fxy{t)\/2. With this choice for P^y, 
we have 

E ">(y-^)^'y/{5(e(:^))+5(C(y))}/(0'd^p. 

x,y£Z'' 

This last term is exactly equal to (r„(^, F{t)), p), which proves that the supre- 
mum in (|7.3p is less or equal than 0. □ 

Notice that the term F) can be written as 

For a given function F, let us define the symmetric part F'^ and antisymmetric 
part F'' of F by 

^x,y — '^{^Px.y + Fy,x), 

^x.y = 2^'^^'^ ~ -Fy,x)- 

We have F) = 83{£_, F'') and r„(^, F) > r„(C, F'^). Therefore, the estimate 
(j7.ip is better for antisymmetric functions For an antisymmetric function F, 
we have 

Let us denote by C^^.^tO^'^) s^*' uniformly continuous, bounded, twice dif- 
ferentiable functions G : R'' x M'' ^ R'' such that G{x, y) = -G{y, x). Notice that in 
particular G{x, x) ^ for any G G Cl^^^{R'^). Let {G'; z = 1, . . . , fc} be a sequence 
of continuous paths G' : [0,T] ^ Cfc2^„t(M'*) and define Fi{x,y) = G^x/n, J//71). 
We will extend the definitions of £„G, CG to antisymmetric functions. For G G 
C,%„t(K'^) we define 

C,,Glix/n) = ^ q(,/^)|Gj(^,^) +Gj(^,^)}, 
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CG{x)^ J q{z){G{x + z,x) + G{x - z,x)}dz. 
With this notation, 

£^i^,Fi) = ^Y1 9ia^))^nGl{x/n). 

Using the second-order Taylor expansion of Gl around {x/n,x/n), we see that 
hm sup \CnGl{x/n) — CGl{x/n)\ = 0. 

These definitions reduce to the previous definitions of £„, C when we take 
G\{x,y) — Gtiu) ~ Gti^)- III the same way, we define 

QnGl{x/n) = \Y1 (li^/n)G{{x + z)/n,x/nf, 



QGix) 



qiz)Gix + z^dz. 



Then, 



r„(^,i^/) = - E 9{m)QnG\{x/n) 



and Q,nG\(xln) is converging uniformly to QG\{x/n) in the sense that 
lim sup \ QnG\{x/n) - QG\{x/n)\ = 0. 

Let the sequence {G\, z = 1, . . . , A;} be fixed. Assume that the trajectories G\ are 
of class . As in Section l473l let tt. be a limit point of the empirical process, and 
let n' be a subsequence such that tt" converges to tt. . Remember that we already 
proved that i:\dx) = u{x, ■)dx for some locally integrable density u[x, t). Using the 
Replacement Lemma in the same way we did it in Section [4.3[ we can prove that 



\dt 



hm sup / {£<:AQ\F\t))~TAif ,F\t))]c 

(l){u{x,t)){ CGI i^) - 2 QGl i^) } dxdt. 
Combining this last line with the energy estimate Ij7.ip . we conclude that 



sup 

l<'i<fe Jo 



E 



sup 

l<i<fc Jo 



<l>{u{x,t)){CGi{x) - 2QGl{x)}dxdt 



By the monotone convergence theorem, we conclude that 



E 



sup 

G Jo 



(l){u{x,t)){LGt{x) -2QGt{x)]dxdt 



< K. 



where the supremum is over the set of continuous paths G : [0,r] C^^^^iK^). 
We conclude that the limit density u{x, t) satisfies 

sup / [ (t){uix,t)){CGtix) -2QGt{x)}dxdt < +00, (7.4) 

G Jo J 
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almost surely with respect to the law of tt.. Let us define (ptix) — 4'{u{x,t)). Let 
II • \\4>.T be the norm defined by 

|2 







||Gm,j, = 2/ I ^t{x)QGt{x)dxdt 

{(j)t{'x) + 4>t{y))h{y ~ x)Gt{y,xfdxdydt 



for any continuous trajectory G : [0,r] — > Cb.ant(IR.'^)- The second definition will 
be more convenient by symmetry reasons. We denote by /^^ rp the Hilbert space 
obained as the closure of such trajectories under || • ||0,t- With this notation, we 
can rewrite (|7.4[) as 

sup|^ j (f){u{x,t))CGt{x)dxdt - ||G|||,t} < 

We recognize in this formula the variational formula for the norm with respect to 
£^ T ant of linear funcional \{G) = 1/2 jj" / (j){u{x, t))CGt{x)dxdt. We conclude 
that A(-) is a continuous functional in C^^Tant- R-issz representation teorem, 
there exists an antisymmetric function Ft{x, y) such that 



(0t {x) + 0t {y))h{y - x)Ft [x, yfdxdt < +oo, and 

KG)^l{\\F + G\\lr~\\F~G\\lr}. 

The arguments above allow to justify the formal computations of calculus of 
variations, from which we obtain that 

/ h[y- X) A t \ : J, f \ dxdydt < +oo. 

Jo J J <f>t[x) + (l)t[y) 

Aside from the factor (j)t{y) + (f>t{x) in the denominator, this estimate is exactly 

the estimate in Theorem 14 .71 At this point we need an extra argument. When g(-) 

satisfies (B), the interaction rate </) is bounded above by ipc and therefore 



£{4>{ut), 4>{ut))dt <24>c / / h{y - x) dxdydt < +oo. 

When g{-) satisfies (FEM), (pc — +oo and another argument is needed. We 
appeal to Hypothesis (C). We point out that when g{-) satisfies (FEM), this 
is the only place where we need to consider an increasing rate function g(-). In 
Appendix IB] we explain how to get rid of Hypothesis (C), but still assuming that 
g{-) is increasing. Under Hypothesis (C), for any function G G Cc(]R'') we have 

limsup7r"(G) < limsup^ ^ p'G{z/n) = p' G{x)dx, 
and in particular u{x, t) < p' for any x € M'', t £ [0, oo). In this case we have 

r £icPiut),^iut))dt<2(t>ip') r ff h{y ~ x) dxdydt < +cx3, 

Jo Jo J J <Pt{x) + 4>t{y) 
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This ends the proof of Theorem 14.71 under both conditions (FEM) or (B). For 
the sake of completeness, now we state a weaker form of Theorem l2.51 which follows 
after checking which hypothesis we have used on each step of the proof of Theorem 
231 



t6 Theorem 7.2. Let g{-) he an interaction rate satisfying sup„ \g{n + 1) — g{n)\ < 
+00 and (FEM). Let p{-) be a transition rate satisfying (P). Let uq : M.'^ [0, oo) 
be a measurable, locally integrable initial profile and let {i^"}n be a sequence of 
probability measures in flzr associated to uq- Let ^" be the zero-range process ^t„ct 
starting from . Assume that there are positive, finite constants p,K such that 

(H) For any n > 0, i?(j^'Vp) < Kn'^. 

Let tt" be the empirical measure associated to Then, the sequence {7r"}„ is 
tight, and any limit point n. o/ {7r"}„ is concentrated on paths of the form u(x, t)dx, 
where u(x, t) is a weak solution of \2.4-^ satisfying (|2.5p . (|2.7p and in place of 

cT 



£{^J(t>{u{-,t)), yj4>{u{-,t)))dt < +O0 (7.5) I ec25.1 







for any T > 0. 



We call these solutions finite entropy solutions, by analogy with the solutions of 
heat equation with finite entropy. This Theorem implies the hydrodynamic limit 
as stated in Theorem 12.51 conditioned on a uniqueness result for weak solutions of 
(123) under ((73)) . 

8. Uniqueness results for the hydrodynamic equation 

8.1. The Unear case. In this section we prove the uniqueness results we need in 
order to establish the hydrodynamic limits of the exclusion process and zero-range 
process with long jumps. We start with the linear case. We learned this proof from 
Luis Silvestre. Let u{x,t) be a weak solution of (|2.ip . By linearity, we can assume 
uo(-) = 0. We will extend u{x,t) to negative values of t by taking u{x,t) = if 
t < 0. For r > 0, we define 9ru{x, t) = u{x, t + r). The function Ot-u is also a weak 
solution of (|2.ip . Since linear, convex combinations of solutions are also solutions 
with the same initial condition, for any integrable function h : [0, a] — > [0,cxd), 

/•a /"OO 

u^{x,t) = i 6T-u{x,t)h{T)dT ~ I u{x,t — s)h{s)ds 







is a weak solution of (12. In particular, if h is the restriction to [0,a] of a 
continuously diffcrentiablc function of compact support contained in (0, a), then 
u^{x,t) is continuously differentiable with respect to time. In the same spirit, let 
h' -.W^ [0, co) be a twice continuously differentiable function of compact support. 
The function 

u''^^\x,t)= [ u''{x-y,t)h'{y)dy 



is also a weak solution of (|2.ip with initial profile uo = 0. But now the function 
u^'^ [x, t) is also twice continuously differentiable in time. Let us suppose that 
u{x,t) is not identically equal to 0. Then, considering h, h' as properly defined 
approximations of the unity, we can assume tha u'^'^ {x, t) is also not identically 
equal to 0. We say that a solution of (|2.ip is classical if the solution is twice dif- 
ferentiable in space and differentiable in time. In that case, (|2.ip is satisfied for 
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each pair {x,t), and not only on a weak sense. We have proved that the existence 
of a weak solution of ()2.1|1 . not identically null and with initial condition = 0, 
imply the existence of a classical solution with the same properties. But for clas- 
sical, bounded solutions, uniqueness is an immediate consequence of the maximum 
principle, as we will see. 

Let u{x,t) be the classical solution constructed above. Fix T > and define 
Ut = sup{u{x,t);x e R'^,t e [0,T]}. If the supremum Ut is attained at some 
point (xojto); we necessarily have Cu(xo,tQ) < 0, with strict inequality if u(-,to) 
is not identically constant. In that case, since (dt — C)u = 0, we conclude that 
dtu{xo,tQ) < and u{xo,-) is decreasing in a neighborhood of tp- This contradicts 
the fact that (a;o,to) is the global maximum of u, unless that = 0. This is 
basically what the maximum principle says. Of course, if u(-,io) is constant or if 
Ut is not attained at any point, we need an extra argument. Let e € (0, 1/2T) be 
fixed, and take {xQ,to) such that u{xQ,tQ) > Ut — £■ Consider the test function 
g{x) = (actually, any smooth function with a strict maximum at x = 

and a fast decay at infinity would be as good as g{x)). Consider A > such that 
f{x) =: g{Xx) satisfies sup^ \£f{x)\ < e and define v{x,t) = 2ef{x — xq) — 2e^i. 
Then, (dt - C)v{x,t) = -2e^ + 2eCf{x - xq) < for any x € R'^, t € [0,r]. 
Therefore, {dt -C){u + v) = {dt - C)u < for any a; G R'^, t e [0,T]. Moreover, 

{u + v){xn,to) > C/t - e + 2e - 2e^to ^ Ut + e - 2e^to > Ut, and 

{u + v)[x,t) <Ut- e^t 

for X far enough from a;o- Now we can conclude that u + u has a global maximum. 
Since {dt — C){u -\- v) < 0, the argument exposed above tells us that the global 
maximum oi u -\- v is attained at t — 0. But u(x, 0) + v{x, 0) = 2ef{x — xq) < 2e, 
so u{x,t) < 2e - v{x,t) for any x G R'^, t e [0,r]. Since v{x, t) <2e- 2eH, we 
conclude that u{x,t) < 2e(l + eT) for any x g R*^, t S [0,T]. Since e is arbitrary, 
we conclude that u{x,t) < for any x,t. Repeating the argument for —u{x,t), we 
conclude that u{x, t) = 0. 

8.2. The nonlinear case. In the nonlinear case, uniqueness of solutions can be 
obtained by an argument due to Oleinik. In order to use Oleinik's argument, we 
first need to prove that solutions have a bounded second moment with respect to 
the reference density p. We start with some elementary lemmas. Here and below, 
c will denote a constant which may change from line to line, and that depends only 
on fixed parameters, like eq, a, d, etc. For simplicity we will assume the ellipticity 
condition ei < (t>'{u) for any u > 0, although the arguments can be carried out for 
any function arising from the interaction rates g{-) considered in this article. 

Lemma 8.1. There is a finite constant c — c{d,a) such that for any b : R"* — > R 
with ||6||oo, ||V5||oo < +00, 

sup / h{y~x){biy)~bix)fdy<c\\b\\l^"\\\/b\\':,. 

Proof. This is an easy application of the mean value theorem. In one hand, {b{y) — 
bix))^ < ||V6||^|y - x\^. In the other hand, {b{y) - b{x))^ < 4||&||^. The first 
estimate is good when \y — x\ is small. The second one is good when \y — x\ is big. 
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After changing to polar coordinates, we have 

//I /•/ A A /"OO 

h{y - x){b{y) - b{x)fdy < / || V6|p^ri^"dr + / ||6|| 

£0 Jo <^o Ji 



eo(2 — a) eoa/" 

where is the area of the unit sphere in R'' and ^ > is arbitrary. Taking 
^ = ||^||oo/|| V6||oo, we obtain the desired bound. □ 

Lemma 8.2 (Poincare inequahty). There exists a positive constant c such that for 
any locally integrable function b ; R*^ --^ R such that £{b, b) < +oo and any M > 0, 

b{xfdx<cM°'£{b,b) + ^(^ j b{x)dxy. 

\x\<M \x\<M 

Proof. In our setting, this inequahty is specially elementary. In fact, it is enough 
to observe that 

{b{x) - b{y)Ydxdy < iM'^+" J J h{y - x){b(y) - b{x)fdxdy, 

\x\<M \x\<M 
\y\<M \y\<M 

and to develop the square in the first integral. □ 

Remember that our objective is to prove that energy solutions of (|2.4|) have 
bounded second moment with respect to the reference density p. The idea is to 
prove that suitable truncations have uniformly bounded second moments and then 
to pass to the limit. Let us define a : R'' ^ R by 

kl<l 

a{x) = {2- \xl 1 < |a;| < 2 
2 < 

and for M > define a^'^ (x) — a(x/M). For a given function u : R'' ^ R, we define 
its truncation u*^ by taking u^\x) — {x)u{x) for any x £ R'^. Let us write 
B{M) = {x eR'^ : \x\ < M}. Notice that u^^ has a support contained in B{2M). 

1 17 I Lemma 8.3. Let m : M'' ^ R 6e locally integrable and such that E{u,u) < +oo. 
Then, there is a constant c independent of u such that for any M > 0, 
is finite and 



£iu^',u'')<c£{u,u) + j^(^ / uix)dx) 



(8.1) I ec25.5 

B{2M) 

Proof. Let us compute £{u^,u'^'^). Since u^' {x) ~ ii x ^ B{2M) and since the 
kernel h{-) is symmetric, we see that 

£{u^\u^^) = JJh{y- x){u{x)a^^{x) - u{y)a^^ {y)f dxdy 

<2 jj h{y~x){u{x)a^Hx)-u{y)a^\y)fdxdy. 

x£B{2M) 



HYDRODYNAMIC LIMIT OF PARTICLE SYSTEMS WITH LONG JUMPS 



35 



Let us write u{x)a'^'' (x) — u[y)a^' {y) = u{x){a^^ {x) — a^' [y)) + a^' [y){u{x) — 
u{y)). Using the elementary inequality (a + b)^ < 2a^ + 2b^, we obtain the bound 



£•(^^*^u*0<4 jj h{y ~ x)u{xf {a^' {y) - a^\x)) dxdy 

xeB{2M) 

ff 2 (^-2) 

+ 4 // h{y-x)a^^{yy{u{y)-u{x)) dxdy. 

xeB{2M) 

Notice that \/a^^{x) = Va{x/M)/M. By Lemma ISTTl the first integral is 
bounded by cM^" Jg^^j^j^u{x)'^dx. Using the fact that a*^ is bounded by 1, the 
second integral in (|8.2p is bounded by £{u, u). Therefore, we conclude that 

£{u'',u^')<j^ J u{xfdx + 4£{u,u). (8.3) 

B(2M) 

The lemma follows using Lemma [5^ to estimate the integral in (|8.3p . □ 



Let u{x,t) be an energy solution of (|2.4p . Let us write Ut{-) — u{-,t) and 
</>*(•) = 0(',i)- Let us denote by (•, •) the inner product in C'^(U.'^). Thanks to 
the energy estimate (|2.6p . relation (|2.7p can be written as 

{ut,dGt)dt- [ £{q^t,Gt)dt+{uo,Go}=0. (8.4) 



When the test function Gt is of the form G{x,t) = g{x, s)ds for a smooth 
function gt-, this relation can be written as 



^ {ut,gt)dt+^£(^J^ <Ptdt,j^ 9tdt)=J^ {uo,gt)dt. (8.5) 

Let us assume for a moment that we can choose gt = (f)t &s a test function. 
In that case the first term in (|8.5p will be comparable to the norm of ut, by 
convexity the second term will be bounded by the energy of Ut and the third term 
will be bounded for suitable initial data uq. Of course, 4>t is not smooth, so we need 
to approximate (j)t in a suitable way. The first thing to do is to truncate 0t. Let 
us define (pf^ = a^^(pt and let {ie',£ > 0} be an approximation of the identity. The 
function (^f^'*^ = * (pt' is a legitimate test function, so all the integrals in (|8.5p 
are well defined. Therefore, we have 

{uu cj^f^dt + ^e* Mt, ^ (Pt'dt) = ^ {i,* uo, cj)t')dt. (8.6) 

Notice that there is a constant c > such that u < c{T-C{u) + 1) for any m > 0. 
By the entropy bound, we have 

T j. 

(t>t{x)dxdt < cT{l + M''). 







B(2M) 
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Therefore, by Lemma 18.31 and the convexity of the energy form, we have 



fdt. 



dt] < 



fdt. 



tdt 



J%rdt,J\fdt)} 



< c 



£{(t>t,(f>t)dt 



(j>t(x)dxdt\ 



B(2M) 



For the term involving the initial condition, we have 



l,UQ,(l)f'')dt < c\\uo\\oo{l + M"). 



Therefore, for M > 1 we have (ut, (ji^ '^)dt < c(l + M'^). By Fatou's lemma, 
{ut,4>f')dt < lim inf e Jq {ut,<l)^'^)dt. We conclude that {ut, (l)f'^)dt < c(l + 
M'^). Notice that > eoUt in B{M). Therefore, 

rT n 

/ u{x, tfdxdt < c(l + M'^) (8.7) 

B(M) 

and in particular u{x,t) belongs to Cf^^lW^ x [0,T]). Let us define ut — ut — p, 
4>t = <t>t — 'Pip)- Relation (|8.5[) can be written as 



{ut,gt)dt + ^£ 



Now we take as a test function 



dt, / gtdt] = 



jM.e 



{uo,gt)dt. 



* (0ta ). We will repeat the computa- 
tions above. Remember that Ti,{u) ~ c{u — pf for u close to p and Ti{u) ^ cu for u 
big. Therefore, there exists a finite constant c such that 7i(u) > cmin{|M — p\,{u — 
p)^} for any u > 0. Therefore, there is a constant c which only depends on T, on 

the entropy bound and on /(wo — pfdx, ||mo|1oo such that Jq {ut, (j)^'^)dt < c. This 
estimate is better than the previous one, due to the introduction of the reference 
density in the definition of (pf 

This time we will use (|8.3I) instead of Poincarc inequality in order to estimate 
the energy form. We have 



£ 



tdt, 



dt] < c 



£{(pt,4>t)dt 



c 



u{xfdxdt. 



\x\<2M 



The energy bound says that we can bound the first integral by a constant c in- 
dependent of e or M. Notice that the second bound now makes sense since is lo- 
cally integrable. Also due to this fact, we have fim^ Jq {ut, (pf'''^)^^ — lo i^t, 4>^}dt. 
Since ut and have the same sign, we conclude that 



ut{xfdxdt < Co 



ut{xfdxdt. 



(8.9) 



B{M) 



B(2M) 



ec27.1 



ec27.2 



ec29 
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We will see that inequality (18. 9p plus the entropy estimate (12. 5p imply that 
/ Ut{xYdxdt < +00. Let us define a„ = Jg^^„^ut{x)'^. By (|8.9p we have 

fln+i > (a„ — co)2""/ci. The idea is that a sequence {a„}„ satisfying this inequality 
is either bounded or has a very fast growth. Let us assume that lim„ a„ = +oo. 
Let no be such that uq > (logj ci + l)/a and a„o > 2co. If n > no, then (a„ — 
co)2""/ci > (a„ - co)2""Vci > 2(a„ - co) = a„ + (a„ - 2co). Therefore, if also 
a„ > 2co, we conclude that a„+i > a„. Inductively we conclude that a„+i > a„ > 
2co for any n > no. Now we start to iterate the inequality in order to get better 
estimates. Since a„ > 2co for n > no, we have Qn+i > 2""co/ci. Therefore, there 
exists a constant (3i > such that a„ > /3i2"" + cq for any n > uq. Inductively, 
let us assume that there is a constant /?p > such that a„ > /3p2P"" + co for any 
n > no. Then, for n > no we have a„+i > /3p2*^^'+^^""/ci and there exists /3„+i > 
such that a„ > /?p+i2(P+^^"" + co for any n > no. 

We have just proved that if a„ converges to +oo as n — > oo, then for any p > 
there exists a constant (3p such that a„ > /?p2P"" for any n large enough. In words, 
the integral of ut{x)'^ over a ball of radius M grows faster than any power of Af. 
But this can not happen. In fact, by (|8.7|) . a„ < c(l + 2"''). Taking p such that 
ap > d and n big enough, we obtain a contradiction. 

We have now proved that any energy solution of (|2.4|) is square integrable with 
respect to the reference density. Let us take two solutions u^, v? of (12. 4|) . Let us 
write = (l){u^{x,t)) for i = 1, 2. By (|8.8p . for any test function gt we have 

^ {ul-ulgt)dt+^£(^J^ (4 - cj)^^)dt, gtdt)^Q. 

But now we know that u^—v? is in £^(R''x [0, T]). Therefore, we can approximate 
u\ — M( by test functions to obtain that 

f I {u\{x) - ul{x)fdxdt + fict^l - cj^Ddt, fic^l - 4>^,)dt) - 0. 

J Q J I/O I' 

Since both terms above are non- negative, we conclude that = u^, and unique- 
ness follows. 

9. The tagged particle problem 

As an application of the results proved above, we obtain in this section the 
scaling limit of a tagged particle in the zero-range process with long jumps. A 
similar result has been obtained for the exclusion process with long jumps [6] and 
for a mean zero, finite-range zero-range process in dimension d = 1 [7]. 

Take a sequence of initial measures {v^} satisfying the hypothesis of Theorem 
12.51 Assume as well that i^"(f(0) > 1) ~ 1 for any n. This condition guarantees 
that the process ^" starts with at least one particle at the origin. Tag one of these 
particles, and follow its evolution, as well as the evolution of the particles as a 
whole. We need to decide in which way the tagged particle interacts with the other 
particles. If we do not want the tagged particle to be different from other particles, 
each time a particle jumps from the site x where the tagged particle is, we decide 
that the tagged particle is the one who jumps with probability 1 {x) . 

We obtain in this way a process (^", X"), where X" is the position of the tagged 
particle at time t. We call A"" the tagged particle process. Notice that due to the 
interaction with other particles, is not a Markovian process. However, it can 
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be shown that X" is a local martingale (it is a martingale for a > 1). We want to 
obtain, under the hypothesis of Theorem 12.51 the scaling limit of the process X". 

Theorem 9.1. Assume that the hypothesis of Theorem \ 2.5\ are fulfilled. Assume 
in addition that d = \, a > \ and that 

(SG) There are positive constants kq, ko such that \g{n + fco) — 9{''^)\ ^ '^o for 
any n € N. 

Then the tagged particle process satisfies 

lim — = Zf 

n — ^oo fi 

in distribution with respect to the J-Skorohod topology o/ 1?([0, cxd), R), where Zt is 
the unique process such that 

exp <^ i0Zt - c\e\" / . „ . rfs } 

[ Jo u{s,Zs) J 

is a martingale for any 9 Cz R. 

The constant c in this theorem is equal to //i(x){e*^ — l}dx. Since we are 
restricted to dimension d = 1, we can assume h{x) = Observe that in 

this theorem a superdiffusive scaling t — > tn" is already embedded in the definition 
of X". The proof of this theorem follows from a careful adaptation of the arguments 
in [5], [7]. All the needed tools we already introduced in the previous sections. We 
give here an sketch of proof, leaving the details to the reader. 

Sketch of proof . Instead of considering the pair we will define an auxiliary 

process Ct by taking C?ix) = Cti^ + ^t") for x 7^ and Ct"(0) = ^("(^D - 1- 
The process is known in the literature as the environment as seen by the tagged 
particle and it was introduced by Kipnis and Varadhan in [8]. Notice that the 
evolution of C" is Markovian. Moreover, the Palm measures i>p{d£,) = P^^(,{0)'^p{d(,) 
are invariant under the evolution of C". The evolution of can be described 
as follows. At each site x ^ 0, a particle leaves x and goes to a; + y at rate 
n°'p{y)g{Q'' [x)) . At site x = 0, a particle leaves to site y at rate n°'p{y)g{C"{0)), 
where g{n) = ng(n + 1) / {n + I) . And finally the whole configuration of particles is 
translated by z at rate n"p(z)6(C"(0)), where b{n) — g{n)/n. The translations of 
the system correspond to jumps of the tagged particle. Let N^''^ be the number of 
translations by z performed by C" up to time t. The number TV^^'" corresponds to 
the number of jumps by z of the tagged particle up to time t. We have the formula 

n n ^-^ 

where Z» Z\ {0}. 

Since A^j^'" is a Poisson process, the process 

Nr- f ny{z)biCnO))ds 
Jo 

is a martingale. Notice that the compensators for zN^'^ and —zN^'^'^ cancel each 
other. After checking the integrability of X", we conclude that X" is a martingale. 
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More convenient will be to consider some exponential martingales associated to X". 
For each S R define 



Alt'" = exp|z6'Xr/n-n" Vp(z)(e'^^/"-l) / 6(C(0))ds|. (9.1) 

Notice that the integrand does not depend on z. The sum over z of the expression 
before the integral is equal to 



n 



which is a Riemann sum for / /i(x){e**^ — l}dx = — c|0|". The proof goes as follows. 
Define the empirical measure tt" by 



nndx)^-J2C{z)S^^„{dx) 

n ^ — ' ' 



n 



As in [6], the sequences {7r"}n, {X."/n}„ are tight. Notice that, aside of a factor 
\/n corresponding to the tagged particle, the empirical measure tt" is the trans- 
lation by X"^ /n of the empirical measure tt" introduced before. We can find a 
subsequence n' such that these two processes converge to some limits. We already 
know that 7r"(c?x) converges to u{x,t)dx. Since the limit u{x,t) is deterministic, 
we have joint convergence of the couple {(tt" , X" /n')}„/ to {{u{x,t)dx, Zt)} for 
some process Zf. We need to characterize the limiting process Zt. Observe that 
tt" (dx) converges to u{x + Zt,t)dx. Notice that the expectation of 6(C(0)) with 
respect to Dp is equal to f3{p) =: 4>{p)/p. In (|9.1[) . we want to replace 6(C"(0)) by 
0(Cr'^"(O))/Cr';"(O), where C^"(0) = (en)-i Elli C(^), since C^^"(0) is a func- 
tion of the empirical measure tt" . Such substitution is known as a local replacement^ 
since it does not involve averaging with respect to a test function. The proof of 
this replacement confines us to dimension d = 1 and a > 1 (see [7, for a detailed 
discussion). In the proof of this local replacement, an spectral gap estimate, uni- 
form on the density, for the process restricted to finite boxes is needed. At this 
point is where (SG) is needed [10]. The only difference with respect to [7] is in 
the proof of the two-blocks estimate, where the corresponding local version of the 
moving particle lemma of Section [5] needs to be invoked. 

After proving this local replacement, we can pass to the limit in (|9.ip to show 
that Alj'" converges to 



Ml = exp \ieZt + c|6i|" / /3(u(s, Z,))ds} 



Notice that (3{p) < k for any p > 0. From this bound we can deduce that the 
sequence of martingales {Ai^'" }„' is uniformly integrable, from where we obtain 
that A4f is a martingale. 

□ 
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Appendix A. The case a = 2 

In this Appendix we explain how to obtain the hydrodynamic hmit for the models 
with long jumps in the case a — 1. For simplicity, we restrict ourselves to dimension 
d = 1; the arguments are the same in any dimension. Let us consider the transition 
rate "pi^y — x)=c/|y — The right scaling for the exclusion process or the zero- 
range process associated to p(-) is v? j logn (still superdifTusive). With this scaling, 
it turns out that the operator £„ defined in ()3.2|) satisfies 

lim sup \Lr,G{xln) - cG"(x/n)| ^ 0, (A.l) 
In fact, let e > be fixed. Remember that 

Let us split the sum in two parts: when \z\ < en and when \z\ > en. When 
\z\ > en, we estimate the term involving G by 4||G'||ooj and we obtain the bound 



nlogn ^ 

|z| >en 



< 



log n ^ e^ log n 

When 1^1 < en, we use a Taylor expansion to write the term involving G as 
G" {xjn^z^ jln} + K^z'^/n'^, where i?" is uniformly bounded. We conclude that 



n log n ^ ' ' L V 7z 

\z\<.en 



< 



< 



G"{x/n)(l-^S^ 



\ogn^^ 



plus a rest bounded by C(G)eV log 71. Since log(7V+l) < Ef=i 1/^ < log(A^+l) + l, 
the right-hand side of the previous inequality goes to when n oo. In order to 
prove the hydrodynamic limits of the processes 77" and ^" , only suitable properties 
for the approximation operators are needed, like the one we just proved. It is 
not difficult to show that £„ fulfills all the properties required. 



Appendix B. Condition (C) and unbounded initial profiles 

In this Section we explain how to get rid of condition (C) when g{-) is non- 
decreasing. Remember that for p < a we have I'p ^ t/^. This is also true for the 
one-site marginals Qp, as proved in Sect ion [2T4l Let us denote by qp^a the coupling 
between qp and qa constructed in Section In other words, qp,a is a probability 
measure in No x No with marginals qp, qa and satisfying qp^cr{x^ < x'^) = 1- For 
simplicity, we will concentrate ourselves in the product initial measures described in 
Section [^31 The idea is to construct a more refined version of the coupling between 
two copies of the zero-range process described in Section 12.41 

Taking a closer look at the coupling (^j,^^) constructed in Section \JM we say 
that ^l{x) is the number of first-class particles at site x at time t. We say that 
£,f{x) — S,l{x) is the number of second-class particles at site x at time t. The idea is 
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that first class particles do not feel the presence of second-class particles, while the 
evolution of second-class particles is modified by the presence of first-class particles 
at the same site. We will consider 4 types of particles: blue, green, red and white 
particles. The dynamics is as follows. We will start with a configuration with no 
green particles, and such that red and white particles do not share a site. Blue 
particles are first-class particles, green particles are second-class particles and red 
and white particles are third-class particles. Each time a red particle jumps over 
a site with at least one white particle, the red particle and one white particle are 
annihilated and a green particle is created. The same happens if a white particle 
jumps over a red particle. 

Let us define this process in a more precise way. We denote by ^f{x), ^^{x), 
^^{x) and ^^{x) respectively the number of blue, green, red and white particles at 
site x at time t. The process (C?, ^f^, ^f", Ct^) is defined in a proper subset of fi^^ 
and it is generated by the operator L = -|- L*^ -|- -|- L, where 

X [fie, {er\ e, r ) - f{e,e,e, r )}, 

X (giei^) + + e'^C^)) - gi^'^ix) + e{x))) X 

x{m'',e,e,{e'rn-fie,e,e,i'^)}, 
if{e,e,e,i'^)= E p{y-xm'^{y)^o)x 

X {g{e{x) + e{x) + ei^)) giei^) + ei^))) x 
X {f{e. e+Sy^e-Sx^e"- Sy) - /r , f, e, n} 

+ E Piy-x)l{eiy)^0)x 

X {g{f{x) + + e'ix)) - 9{e{x) + fix))) X 

X {f{e,e + Sy, e - Sy, r - Sx) - f{e,e, e, r )}• 

The first line govcrncs the evolution of blue particles and it is independent of the 
other particles. The second line governes the evolution of green particles, and it is 
independent of the evolution of red and white particles, except for the creation of 
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new green particles, which corresponds to the last two lines. The third and fourth 
lines correspond to the evolution of white and red particles when they do not meet. 

Let Mo(-) be an unbounded profile, satisfying the hypothesis of Theorem 12.51 
Let us recall the definition of /i": the measure /i" is the product measure in flzr 
of marginals given by ^"(^(x) — I) — qu^(x){l)- Fix M > 0. For each n, let us 
consider the product measure p." in defined by 

where u^'^ = minjup (x), M} and C/*^'" — maxj-Ug (x), M}. Now we define the 
initial measure /i" in Slf^ by taking 



0, if uo(a;/n) < Af 

e{x)-e{x), iiuo{x/n)>M, 



i^{x) 



^2(a;) if uo(a;/n) < M 
0, if uo{x/n) > M, 



After a careful checking, we see that evolves like a zero-range process with 
initial distribution /i"'^^ associated to the profile min{uo(a;). A/}, + + Cf" 
evolves like a zero-range process with initial distribution ^" and that + + 
evolves like a zero-range process with initial distribution I'm- Therefore, Theorem 
12.51 applies for the rescaled process . Sending M to infinity and using the 
monotone convergence theorem, we see that the rescaled process ^(„q with initial 
condition ju" also satisfies the conclusion of Theorem 12.51 



Appendix C. Condition (H) and bounded profiles 

In this section we explain how to get rid of condition (H) when g{-) is non- 
decreasing and the initial profile is bounded. Let mo(') be a bounded profile and 
let {/i"}n be the sequence of measures defined in Section [2T4l Assume that there 
are constants < po < Pi such that po < uq{x) < pi for any x € R*^. We will use a 
coupling similar to the one used in the previous section. The idea is the following. 
Take A/ > and define 

Uq {x) 



uo{x), 


\x\ 


< M 


Pi, 


\x\ 


> M, 


uo{x), 


\x\ 


< M 




\x\ 


> M, 



■ {x) 

We have Wq^' (x) < uq{x) < Uq^'^{x) for any x G M''. We will have three classes 
of particles: blue, green and red particles. The initial configurations are as follows. 
For any three numbers Po ^ P !i Pi it is possible to construct a measure qp„^p^p^ in 
N X N X N with marginals Qp,,, Qp, Qp^ and such that q{x^ < x^ < x^) — 1. For each 
n, let us consider the product measures //" in 17^^ defined by 

~^-{^^{x)^lB,e{x) = lG,e{x)=lR)^ 

= 1u'^'-{x/n).Uo{x/n),u^''-+(x/n)i^BjB + h, Ib+Ig+ Ir)- 
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The process {S,^ ,£,^, with initial distribution /i" is generated by the operator L 
defined in the previous section, projected on the set of configurations without white 
particles. As we pointed out before, evolves like a zero-range process with initial 
distribution associated to the profile Uq''~ , +Ct^ evolves like a zero-range process 
with initial distribution uq and + ^f' + evolves like a zero-range process with 
initial distribution u^''^. Both rescaled processes and £,^a +£,^c satisfy 

the hypothesis of Theorem 12.51 Tightness of the empirical measure associated 
to + follows by comparison, and the limiting measures TTt{dx) — u{x,t)dx 
(possibly random) are bounded between u'^'~ {x,t)dx and u^'~^{x,t)dx, which are 
the solutions of the hydrodynamic equation (|2.4p with initial conditions Uq '~ , Uq 
Sending M to oo we obtain that the rescaled process -I- also satisfies the 
conclusion of Theorem 13.11 When uq is not bounded below, a similar, but more 
sophisticated coupling can be constructed. 
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